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A.  INTRODUCTION 


This  is  the  Final  Report  for  Office  of  Naval  Research  Contract 
Number  N00014-77-C-0586.  This  report  describes  theoretical  studies 
related  to  man-made  disturbances  of  the  upper  atmosphere  which  can 
affect  VLF  wave  propagation.  The  purpose  of  our  research  was  to 
estimate  the  effect  that  artificially  created  plasma  holes  could  have 
on  whistler  made  signals.,  propagating  through  the  ionosphere.  In  the 
course  of  our  research,  we  K'ave  developed  models  which  describe  the  flow 
of  reactive  vapor  which  can  be  injected  into  the  plasma  layers.  Our 
research  has  considered  the  influence  of  the  modified  ionosphere  on 
propagation  of  HF  and  VHF  as  well  as  VLF  radiowaves.  During  the  final 
period  of  research,  we  have  found  evidence  of  man-made  perturbation  of 
the  earth-ionosphere  waveguide. 
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B.  NEUTRAL  GAS  DIFFUSION  MODEL 


An  overview  at  the  status  of  our  ionospheric  depletion  modeling 
indicated  that  our  simulations  of  neutral  gas  expansion  needed  the  most 
improvement.  Consequently,  we  have  developed  a  model  of  three-dimensional 
expansion  of  neutral  gases  in  a  nonuniform,  chemically  reactive  medium. 
This  model  is  described  in  the  following  paper  which  was  published  in  the 
Journal  of  Geophysical  Research  and  which  was  presented  at  the  1978  Summer 


Computer  Simulation  Conference. 
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Three-Dimensional,  Time-Dependent  Modeling  of  Neuiral  Gas  Diffusion 
in  a  Nonuniform,  Chemically  Reactive  Atmosphere 

Pali  A.  Bi:kn!iakdi 


RiiJtii'dcii.i  iMhoratory.  Stanford  L  niremiv.  StatifurJ,  California  Ml' 


A  time-varying  model  .if  peulval  g.is  expansion  m  a  nonumloi m  environment  is  developed.  The  model 
includes  diffusion  in  a  multicomponent  atmosphere,  chemical  reactions  between  the  diffusing  paces  and 
the  atmosphere,  thermal  diffusion  elTects,  and  transport  due  to  altitude-dependent  winds  The  three- 
dimensional  dill'usion  equation  governing  the  neutral  gas  floss  is  solsed  numerically  using  Fourier 
transform  and  It  title  difference  techniques,  examples  of  1 1,.  OH.  and  CO,  dill'usion  illustrate  the  effects  of 
chemical  reactions  anil  ssind  shears  on  the  neutral  expansion.  The  model  may  he  applied  to  chemical 
releases  sshich  produce  ionospheric  depletions  or  luminescent  trails. 


1.  Intromit TION 

Passive  experiments,  such  as  tracing  neutral  ssind  patterns, 
or  active  experiments,  such  as  chemical  depletion  of  the  F 
layer,  involve  the  release  of  neutral  gases  into  Ihe  upper  atmo¬ 
sphere.  Diffusive  expansion  governs  the  floss  of  these  gases. 

Previous  theoretical  treatments  have  considered  various  ap¬ 
proximations  to  diffuse  e  gas  expansions.  An  analytic  solution 
for  diffusion  of  a  point  release  in  a  lossless,  exponential  atmo¬ 
sphere  ssas  presented  by  I  u  and Klein  [  I9n4],  This  solution  svas 
generalized  to  continuous  releases  and  extended  for  height- 
independent  chemical  losses  by  Bernhardt  ri  ai.  [1975]  and 
Bernhardt  [  197t«j.  Forhes  and  Meo-nllo  [1976]  and  Mendillo  and 
Forhes  [  I97S]  have  produced  an  an. .lytic  solution  for  dill'usion 
in  an  atmosphere  w  illi  a  chemical  loss  '.hat  varies  exponentially 
with  altitude. 

A  numerical  model,  as  opposed  to  an  analytic  model,  of 
neutral  gas  expansion  in  a  cylindrical  geometry  has  been  de¬ 
scribed  by  Hairier  and  Carr  [1971],  This  model  described  the 
solution  to  the  diffusion  equation  in  a  two-dimensional  (radial 
and  axial)  coordinate  system  using  a  finite  difference  tech¬ 
nique.  This  model  gives  more  realistic  results  than  the  analytic 
models  because  it  does  not  require  a  simple  analytic  descrip¬ 
tion  of  the  atmospheric  constituents.  However,  the  two  dimen¬ 
sionality  of  the  model  limits  the  range  of  its  applicability. 

None  of  Ihe  models  mentioned  thus  far  have  included  (I  I 
the  effects  of  an  altitude-dependent  neutral  wind.  (2)  the  diffu¬ 
sion  of  neutral  species  produced  by  chemical  reactions  be¬ 
tween  the  released  gas  and  the  background  atmosphere,  (.1)  the 
expansion  of  several  gases  (which  may  react  with  each  other  as 
well  as  with  the  atmosphere).  (4)  the  effects  of  thermal  diffu¬ 
sion.  and  (5)  the  effects  of  a  multispecies  background  atmo¬ 
sphere  This  paper  describes  a  numerical  model  of  dill'usion 
which  contains  all  of  the  features  mentioned  above. 

I  he  numerical  model,  however,  has  its  limitations.  Latitudi¬ 
nal  and  longiludin.il  variations  in  the  background  atmosphere 
ai’d  neutral  wind  velocity  are  not  included  Also,  any  depar¬ 
tures  from  pure  diffusion  are  neglected.  I  or  example.  Sehttnk 
[ITS]  has  shown  that  nonlinear  acceleration  and  viscous 
sties'  may  become  significant  lor  releases  above  500- km  ulii- 
ludc.  A  tnimerie.il  model  including  these  effects  is  currently 
being  developed 

2.  l  lIKI  I  -Divtl  NSIIINXt  Dim  SION  I. CA  MIONS 

I  he  equations  describing  ihe  dill'usion  of  several  gases  in  a 
multicomponent  atmosphere  are  developed.  I  ourier  trunsfor- 
(  npvrn'tii  c  IT'J  by  the  American  Geophysical  Union 


mation  is  used  to  remov  e  the  spatial  derivatives  in  the  horizon¬ 
tal  plane.  An  exponential  transformation  is  used  along  the 
vertical  axis.  In  the  formulation  the  subscripts  /'  and  k  refer  to 
the  r'th  injected  gas  species  and  the  Ath  atmospheric  species, 
respectively.  The  coordinates  v.  y,  and  :  represent  the  east, 
north,  and  vertical  directions,  respectively. 

Consider  M  neutral  species  diffusing  through  a  background 
atmosphere  composed  of  A'  gas  species.  The  diffusion  velocity 
for  neutral  species  i  is  written  as 


w,  =  w  -  l), 


-LV„,  +  A.q<'  + 

n,  //, 


vrj 


(I) 


i  =  1.2. 


.  M 


[I'ineenti  and  Kruger,  1965;  Banks  and  Kockarts.  1973].  where 
the  tth  neutral  gas  has  velocity  vv„  scale  height  II,  =  (A  T) 
(tn,g)  (A  is  Boltzmann's  constant),  mass  nt„  diffusion  coeffi¬ 
cient  />,.  concentration  n„  and  thermal  diffusion  factor  The 
neutral  atmosphere  is  characterized  by  a  temperature  7,  grav¬ 
ity  acceleration  gt i,.  and  velocity  w.  1  he  unit  vector  in  the  r 
direction  is  a,.  The  diffusion  coefficient  is  found  from 


1 

'ri, 


(2) 


where  l>,„  is  the  mutual  diffusion  coefficient  between  the  r'th 
and  the  Ath  neutral  gases.  The  effect  of  the  injected  gases  on 
the  diffusion  coefficient  is  neglected.  Fquation  (I)  is  valid 
only  when  the  diffusing  gas  is  a  minor  constituent. 

The  continuity  equation  for  the  r'th  neutral  species  is 


4  V •(„,«•,)  =  s,  (31 

where  r  is  tune.  /’,  and  /.,  are  chemical  production  and  loss 
terms,  respectively,  and  S,  is  the  sou.ve  function  of  gas  ».  The 
background  atmosphere  is  assumed  to  be  horizontally  strati¬ 
fied  Consequently,  the  neutral  wind  velocity  (w).  the  diffusion 
coefficient  (/>,).  and  the  neutral  temperature  (7)  vary  only 
along  the  altitude  coordinate  ( - ).  Under  these  conditions  the 
diffusion  equation  is  given  by 
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f  tAtt  *  Ti:  *  <>  "] +  '•  - c  -  i. 

(4) 

The  v  ami y  derivatives  are  removed  by  Fourier  transforma¬ 
tion.  The  two-dimensional  I'ourier  transform  of  n,  is 

A,  -  J  J  n,(x.  y.  I)  exp  +  />>')]  dx  dy 

In  the  r  direction  an  exponential  transformation  is  applied. 
The  exponential  transformation  is  given  by 

?=-//,  exp 

where  :0  is  the  release  altitude  and  II,  is  a  constant.  The 
combination  of  the  two-dimensional  Fourier  transform  and 
exponential  transform  will  be  referred  to  as  the  Fourier-ex¬ 
ponential  transform  (FET).  With  this  transformation.  (4)  be¬ 
comes 
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A,  -  A*,  +  jXi,  C,  -  GHi  +  jG,,  (8) 

Substituting  (7)  and  (8)  into  (5)  and  equating  the  real  and 
imaginary  parts  yield 
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m  -  1 

plus  a  similar  equation  with  XH,  replaced  with  Nu,  Af/i  replaced 
with  -A*,,  and  C,H,  replaced  with  G„. 

The  coupled,  second-order  differential  equations  governing 
the  expansion  of  several  gases  in  a  nonuniform,  chemically 
reactive  environment  are  summarized  in  terms  of  the  matrix 
differential  equation 
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(10 b) 

A,  O' 

A , 

A  = 

(10c) 

where  <5\.  <C„  and  S,  are  the  FET's  of  Pt,  L,.  and  5,.  respec¬ 
tively. 

To  gain  insight  into  why  the  exponential  transform  is  used, 
consider  the  following.  The  diffusion  coefficient  D,  increases 
with  altitude  The  variable  l2  decreases  with  altitude.  Con¬ 
sequently.  the  factor  (0, $*)///,’.  which  multiples  8,S \/8l2, 
may  be  made  relatively  constant  in  altitude  with  an  appropri¬ 
ate  selection  of  II,. 

It  is  convenient  to  partition  the  production,  loss,  and  source 
functions  as  follows: 

P,  -  +  S,  =  £  y.  r.  <)  -  g,(x.  y.  r)  (6) 

m  »  I 

where  /■’,«,  is  an  altitude-dependent  function  which  represents 
the  chemical  coupling  of  species  n„  to  the  production  or  loss  of 
species  n,.  Fourier-exponential  transformation  yields 

(V,  •  if,  +  A,  -  j£  F,m$„  -  CAj,.j,.l)  (7) 

m  ■  I 

where  (/,  is  the  III  of  fi(. 

I  he  complex  quantities  A,  and  (t,  are  written  in  terms  of  real 
and  imaginary  parts; 
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Fig.  I.  Diagram  of  the  numerical  solution  of  the  diffusion  equation  in  three  dimensions.  FF7  stands  for  Fourier- 

exponential  transform  (see  text). 


The  individual  elements  in  the  2\f  X  2 M  matrices  are 
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arid  is  defined  by  (7). 

A  few  words  about  these  equations  are  in  order  The  real 
and  imaginary  parts  of  A  ,  are  coupled  only  by  the  parameter  C 


in  the  B  matrix.  The  parameter  results  from  nonzero  horizon¬ 
tal  winds.  !f  these  winds  are  altitude  independent,  the  transfor¬ 
mation 

S,  =  St  exp  >2irt fxwx  +  /vuy ) z/zj 

G,  =  G,  exp  fjlir (fx*x  +  />■,)<//! 

results  in  a  matrix  with  C  =  0.  This  transformation  is  equiva¬ 
lent  to  following  the  gas  diffusion  in  a  reference  frame  which 
moves  with  the  wind. 

The  equations  presented  in  (10)  have  several  properties 
which  aid  in  their  solution.  The  partial  differential  equations 
are  second  order  in  {  and  first  order  in  l.  The  coefficient 
matrices  A  and  £  are  diagonal  and  do  not  depend  on  S'. 
Coupling  between  the  dements  in  S  occurs  via  multiplication 
with  the  B  matrix  and  via  the  G  matrix.  Since  6  is  independent 
of  S.  the  only  nonlinearity  in  (10)  occurs  through  G. 

In  the  solution  of  (10),  however,  the  nonlinearity  provided 
by  G  may  he  neglected.  The  relative  importance  of  the  nonlin¬ 
earity  varies  according  to  the  production  and  loss  process 
represented  by  O'.  Complicated  forms  of  G  arise  from  chemical 
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OH  diffusion  coefficients  are  found  from  O  diffusion  coefficients  bv  =  0.9141(1  »  (I  f  in  m  tUcKrn-ie  el  al 

1973).  and  (»)  10'.  Read  8.31(17)  as  8.31  x  10".  ' 


reactions  between  two  different  diffusing  gases.  This  chemical 
coupling  is  generally  weak  for  the  gases  of  interest.  Con- 
scquently,  G  is  assumed  to  be  a  function  of  A’  at  the  beginnine 
of  each  step  in  the  solution.  With  this  approximation.  (10)  is  a 
set  of  linear,  coupled  partial  differential  equations.  Before 
these  equations  can  be  solved,  the  initial  conditions  and  the 
boundary  conditions  must  be  specified. 


3.  Gas  Rtt  txsr  Solkc  t  Functions 
and  Initiai.  Conditions 


The  injection  of  a  neutral  gas  may  be  incorporated  into  the 
simulation  either  as  ( I )  an  initial  condition  for  point  sources  or 
short  extended  sources  or  as  (2)  a  source  term  S,  (see  (3))  for 
long-duration  extended  sources.  Point  sources  are  considered 
first  During  early  times  (time  of  release  less  than  one  collision 
time)  a  tenuous  gas  cloud  will  undergo  free  molecular  expan¬ 
sion  The  neutral  gas  concentration  in  space  and  time  for  this 
type  of  expansion  is  given  by 


n(.v.  y.  r,  /) 


/  -(.v*  +  r*  +  ;*)Y 

Aexp  — (^ji- —  ’J 


[Bicnkowski.  1964],  where  ,V„  is  the  number  of  molecules  re¬ 
leased.  (•„  -  (2K  3  is  the  speed  of  sound  inside  the  gas  at 

the  time  of  release,  7„  is  the  initial  gas  temperature,  m  is  the 
molecular  mass,  and  r  is  the  time  since  release. 

At  later  times  (typically  10  collision  times  or  greater),  diffu¬ 
sive  expansion  applies.  The  solution  for  diffusive  expansion  of 
a  point  release  in  a  uniform  atmosphere  is 


1  '  (4rr/70/)’ 1 

where  />„  is  the  diffusion  constant  (e  g  .  Josl.  I%0], 

Both  the  initial,  free-molecular  expansion  and  the  later  dif¬ 
fusive  expansion  have  the  form 


n„(  v.  y.  :) 


-V.  exp  [-  (y!_T/JL -Y'i-ff'Jl 

(2irrr’)’  * 


(13) 


where  a  is  the  ‘radios'  of  the  expanding  gas  at  the  start  of  the 
simulation.  Consequently .  this  expression  is  used  for  the  initial 
distribution  of  a  point  release  in  the  numerical  simulation  All 
simulations  of  points  releases  are  started  three  atmospheric 
collision  times  after  release. 

A  line  source  is  considered  next  To  model  a  source  moving 
in  a  straight  line,  one  can  break  the  release  into  an  infinite 
sequence  of  point  releases,  this  is  done  by  Bernhardt  ( 1976] 
with  the  result 


n,( x .  i .  l) 
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where  f  0  is  the  flow  rate  of  the  source:  is  the  duration  of  the 
source:  r„  is  the  time  that  the  gas  takes  to  thcrmalire  with  the 
background  atmosphere:  D,  is  the  injected  gas  diffusion  coeffi¬ 
cient:  l'..  f\.  and  are  the  velocity  components  of  the 
moving  source;  and  Xc.  Yc.  and  Zc  are  the  centers  of  the  line 
source.  Each  segment  of  the  line  release  is  delayed  by  time  r„, 
which  is  taken  to  be  three  atmospheric  collision  times. 

The  long-duration,  extended  source  requires  the  continual 
injection  of  material  during  the  time-varying  simulation.  To 
accomplish  this,  the  source  term  .S’,  in  the  continuity  equation 
(3)  is  used.  The  source  term  for  a  line  source  has  the  form 


5,(.v.  v.  r.  /)  -  “Tp  exp  j-  [A  -  V, ^ y  -  t'j  -  Xc 

-  -  r,(i- -  -  z,]’}/2"' 


where  the  variables  have  already  been  define  .  in  (12)  and  (13). 
It  is  computationally  easier  to  include  the  gas  source  as  an 
initial  condition  rather  than  as  a  source  term  in  the  continuity 
equation. 


■4.  Bot  NDVR)  Conditions 

The  solution  of  any  second-order  partial  differential  equa¬ 
tion  requires  boundary  conditions.  The  boundary  conditions 
at  the  horizontal  border  of  the  solution  space  are  implicitly 
specified  by  using  the  two-dimensional  finite  Fourier  trnns- 


1  AIM  1  2  Relative  Importance  of  Thermal  Diffusion 

(I  -t  o,l  0  In  7  r:. 


Altitude,  km 

I/f/..  m  1 

m  * 

150 

3.5  x  10  " 

7  X  10  ■ 

225 

2  7  a  |()  • 

7  X  |0  ’ 

300 

2b  X  III  * 

8  X  10  1 

375 

2.5  x  10  * 

8  x  to  ’ 

450 

2.4  x  It)  » 

8  X  to  1 

'WdtL'r  •*  \ 
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TABLt.  -V  Ncuu.il  (i.is  RcIcum*  Simulations 


Simulation 

Species 

Amount 

Released. 

molecules 

Altitude. 

km 

T >  pc  of 
Release 

Wind 

Velocity. 

m/s 

f-  ipure 
Reference 

i 

-o 
=  ° 

T  127 1 
■V  |27J 

300 

point  source 

0 

3-6 

2 

Hj 

-V(27| 

250 

line  source 
r,  -  8  k  in 's 
tu  -  5  s 

0 

7 

} 

CO. 

3.[27] 

160 

point  source 

Given  in 
Figure  8 

9 

4 

CO, 

3  [27) 

160 

line  source 
vM  -  1  km/s 
i'i  -  1km  s 
^  30 s 

Given  in 
Figure  8 

10 

Note  th.it  M  -  10'.  Read  3.[27]  us  7  X  10”. 


form  This  transform  assumes  that  the  solution  is  periodic  in 
tw o  dimensions  with  periods  equal  to  the  horizontal  dimen¬ 
sions  of  the  solution  space.  Care  must  be  taken  to  insure  that 
the  solution  does  not  expand  into  these  borders. 

The  upper  and  lower  boundary  conditions  to  the  second- 
order  partial  differential  equation  in  r  and  t  are  specified  by  the 
equation 

A,  4  H0.\  -  <70  °  (15) 

A  simple  condition  would  be  to  force  the  densities  to  zero  by 
letting  =  0. 


0 

I 


and  C o  -■  0.  This  condition  requires  an  infinite  sink  of  particles 
at  the  boundaries. 

Another  possible  condition  requires  that  the  flux  of  parti¬ 
cles.  po  to  zero  at  the  boundaries.  From  (I  i. 


-  n, ic.  -  l),  -y •  +  (I  +  a,)  -  D,~—  ~  0 

//,  r:  i2 


With  this  assumption  (and  the  exponential  transform)  the 
matrix  elements  in  (15)  become 


A°- ir, D' 


Do,  =  «,  4  /), 


~ff,+  yf;(1  +" 


Co,  -  o 


This  boundary  condition  does  not  allow  any  particles  to  lease 
the  solution  space. 

The  boundary  conditions  described  above  require  infinite 
velocity  and  zero  velocity,  respectively,  at  the  top  and  the 
bottom  of  the  solution  space.  A  condition  between  these  two 
extremes  is  used  in  the  present  simulation. 

The  diffusion  velocity  at  the  borders  is  specified  by 

H,  =  (r  -  (16) 

where  r„  is  the  altitude  of  release  at  time  t  -  0.  Equation  (16) 
implies  that  the  average  velocity  that  a  particle  has  after  mov¬ 
ing  a  distance  c  -  r0  in  time  /  is  half  the  velocity  that  it  would 
have  if  it  encountered  no  collisions.  This  is  indeed  the  gas 
velocity  for  diffusion  in  a  uniform  atmosphere  [ Bernhard i. 
19761.  ' 

Substitution  of  (16)  into  (I)  yields 


/ 

ri.yv.  ■ 


..  iLfii j 

2 1 


"■  [zz: 


+  (■+«,) 


P  In  T 
P: 


l  _  Pn 

>~  D,  —  =  0 


jco  o  cube  centred  at  300.0  altitude 
Tiyf  AFTER  B£i(aSC:  31. t  •« 

*,  CO^CEnTRA*  at  S  jRFACC:  3.0(l|  cm  1 


300  0  kw  CUBE  CfNTfPEO  AT  300  0  km  AlTlTuOE 
Time  AFTER  RELEASE:  31  9  ik 
0m  C0NCENTPAT  ON  AT  SURfACE:  t.9(J|  cm*1 


300.0  km  CUBE  CCNTtREO  AT  JCO.O  km  AUlFUOE 
TIME  AFTER  RELEASE  3i.t  ik 
CO,  CONCENTRATION  AT  SURFACE:  8  3(8)  cm’* 


I  \c  '  StfTUiSt.incoiiN  expansion  of  3  >  HP7  molecules  of  CO,  and  3  y  10”  molecules  of  H,.  released  at  300- km  altitude. 
The  reaction  between  O  and  II,  produces  the  Oil  radical  At  early  times  the  expansion  is  nearly  spherical 


hi  KSH  \kd  i  Mui>n  ini.  n»  Nn  i  k  \  i  C  j  \s  Dun  sms 


300  0  km  Cu9C  C£nT£R£0  AT  300.0  km  AUUuDt  300  0  km  Cu9E  CCnT£R£0  AT  3  DC  0  km  AllrToCf 

T.m£  *rr(4  ff(l£ASC  22 9.3  T.u£  tfl[#  RELEASE:  229  3  <»e 

h j  CO*«C£N’flAT.CN  AT  SuRFACC.  2.2(51  cm'*  CO,  CONCENTRATION  AT  SURFACE:  7.»£ 5 j  cm  1 

Kiji  4  l  .r.c  lime  expansion  of  J  V  10*’  molecules  ofC  0;  jn«J  1!,  .it  300- km  altitude.  1  he  nonsphenc.il  expansion  is  due  Co 

the  altitude  variation  in  the  b;..k-’ round  atmosphere. 


With  these  bour.il.ir>  conditions  the  solution  to  (II)  is  com- 
pitted  mimeucallv 


5.  Sol  fllON  Of  lilt  TlIRH-DlMl  SSttINAI 
Dim  sion  bgiMioss 

1  he  procedure  for  the  solution  of  (10)  is  displayed  in  figure 
I  I  ach  diffusing  p.is  concentration  n,  and  function  y.  is  stored 
in  a  memory  block  w  ith  dimensions  corresponding  to  the  x.  t . 
and  :  ct'ordm.ite  axes.  I  he  gas  concentration  block  is  filled 
with  either  initi.il  values  (at  the  start  of  a  simulation)  or  values 
from  the  solution  at  a  previous  time  (/,)•  The  source  function 
.V,  is  determined  from  external  inputs.  If  the  source  function 
has  an  analytic  f  ourier  transform,  the  numerical  fast  f  ourier 
transform  (II  I  )  will  not  be  nccess.iiy  l  or  instantaneous 
releases  the  source  function  is  expressed  as  an  initial  condition, 
and  .S’,  is  neglected  (see  section  3) 

Horizontal  slabs  in  these  blocks  arc  l  onricr-exponcntial- 
tr.insformeil  and  repacked  inlvi  memciry  blocks  corresponding 
to  f , .  /.,  and  (  coordinate  space  I  he  elements  of  the  produc¬ 
tion  loss-source  function  ((»*,.  (>/,)  arise  from  chemical  reac¬ 
tions  between  the  Indiv  idu.il  diffusing  space  with  an  additional 
term  vine  to  external  sources  I  he  i.  /(.  and  I.  matrices  arc 
determined  from  the  background  atmospheric  parameters 
Vertical  columns  in  are  retrieved  from  data  memory 
blocks  at  fixed  values  of  f,  and  These  are  used,  alone  with 
the  partial  differential  equation  given  in  ( 10).  to  find  values  of 
,V„,  and  ,\„  at  later  times.  The  coupled,  vector  set  of  linear, 
second-order  (parjhotic)  partial  differential  equations  are 
solved  hv  eerier. ili/ing  the  method  ol  Crunk  and  \ut>Ut>n 
|I'U7] 


The  new  solution  vector  is  placed  in  an  output  memory 
block  in  the  form  of  a  •column'  at  /,  and  /,.  The  process  is 
continued  until  the  updated  columns  have  been  determined  for 
all  spatial  frequencies. 

The  inverse  IT  T  converts  each  slab  in  />,  and  £  space 
into  one  at  the  x,  y.  and  e  spatial  coordinates.  T  he  output 
memory  block  contains  the  neutral  gas  densities  at  time  (r,  + 
At)  The  solution  process  is  repeated  using  lime  increments  At 
which  arc  proportional  to  I  'r. 

In  the  above  description  it  has  been  conceptually  convenient 
to  consider  four  separate  memory  blocks:  (1)  .x.  y.  :  space  at 
time  (21  f,.  />.  £  space  at  time  t,.  (3)  />.  £  space  at  time  t, 

-  At.  and  (4)  x.  y.  ;  space  at  time  f,  +  At.  Actually .  these  four 
blocks  should  use  the  same  memory  locations  to  reduce  data 
storage  requirements. 

Because  of  the  uncertainties  at  the  boundaries  (see  section  4) 
an  adaptive  grid  is  adopted  As  the  concentration  oT  the  neu¬ 
tral  gases  becomes  significant  near  the  edges  of  the  solution 
grid,  this  space  is  increased  by  redefining  the  spacing  of  the 


i  to  too  tooo 

TIME  AFTER  RELEASE  (J#t) 


1  ie  s.  Temporal  variation  in  the  maximum,  central  concentra¬ 
tion  of  the  expanding  CO.  and  II,  gas  clouds  T  he  point  of  maximum 
concent-ation  drops  in  altitude  with  tunc  Ihe  concentration  at  this 
potm  i-  greater  than  the  gas  concentration  for  expansion  in  a  uniform 
atmosphere. 


Hi  ksh  \  hi»i  Muiiiiim,  mi  Nu  mi  li\s  Hiiiimmn 


7  mi 


IK'  imiUi.il  diffusion  an ilictcnts  fm  the  diffusing  f.r.c% 
(C  O  I'  ..ml  Ol I )  arc  fou (ul  from  (he  ruruml.i 

/),,  A./I'-n. 


I  i t  I  er  por.i!  v an. 1110:1  1::  the  rn.'.summi  (Ml  ivircc:ilr.itu>:i 
1  he  nuli.il  i:-e  .  cun.cntrjlk'ii  is  ..iu-,.‘J  Hs  reaction  beUsccn  atomic 
os;,  gen  ,i:ui  the  released  I  h  I  he  t.ill  1:1 1  lie  eoiieenlr.il  k'11  is  due  111.mil> 

to  sKKsiea 


sshcre  /  is  the  11eal1.il  tcmpci  al'ire,  n„  is  the  conccit'.iatkin  of 
the  Atii  .ilituisphei le  species.  and  I,,  .mil  .S,,  are  given  in  Table 
I  I  he  diffusion  coefficient  lor  OH  is  found  hj  s.ahng  lire 
diffusion  coefficient  for  O.  according  to  their  lel.iti-.e  masses 
(see  toolnole  to  lahle  I).  The  iHiiiierie.il  sallies  fo ;  ali  diff’ii- 
sum  coefficients  are  hased  on  laboralors  nieasureii'enis  tabu- 
laled  hs  Mci'on  di'.J  Marrero  [1970]  ami  Hank  1  nnJ  A  i r! ' 
|1971j  Ion-neutral  collisions  have  less  than  a  0  |‘  effect  an 
the  diiTiisunt  and  are  therefore  neglected. 

I  he  sallies  of  the  therm.il  diffusion  factors  are  risen  as 
lolhnss  «i ,, ,  -  0.25,  n,,H  0.0.  and  it,  -r  0. 1 .  I  he-e  factors 
ssere  tficorcticalh  derisecl  hs  (  koj’inon  and  ( ’tm/tng  [I970J. 

In  regions  of  large  teinperaluie  gradients  the  effects  of  ther- 
ni.il  di.'Tiision  become  important  A  comparison  of  the  factors, 
sshi.ii  multiple  n.  and  i  n,,'iz  in  (4),  shosss  that  thermal  diffu¬ 
sion,  should  he  considered  if 


(I 


fjn  r 

(■: 


end  Thus  the  solution  grid  follosss  the  outer  edges  ol  the 
doftisis e  esp.i nsion.  Csiiitainine  the  gases  ss  ithin  its  boundaries 

b.  Soil  1  ation  01  C  O...  H..  ani>  Oil 
Dll  1  c  sis  1  Laivsnsions 

I  he  procedure  outlined  in  the  pres  ions  sections  is  used  to 
model  carbon  dioxide  and  or  molecular  hsdrogen  expansions 
C  O.  ami  II.  ssere  chosen  because  the  large  difference  hetsseen 
their  nudectil.tr  masses  produces  a  noticeable  clillcten.ee  in 
their  rates  of  expansion.  Hsdrogen  is  of  special  interest  be¬ 
cause  it  reacts  ssith  atomic  oxsgcn  in  the  atmosphere,  piodnc- 
mc  the  hsdroxsl  radical.  All  three  substances  (CO...  II..  and 
OH)  cause  plasma  depletiims  sshen  injectcsl  into  the  A,  laser 
ionosphere  [Hcrnkortlt,  197b] 

1  lie  neutral  atmosphere  used  in  the  simulation  is  the  empiri¬ 
cal  global  model  obtained  from  mass  spectrometer-incoherent 
scatter  (MSISi  data  hs  UcJtr.  ,7  ol  (1977,;.  h |  l  igtire  7  illus- 
tr.ites  the  composition  and  tenipeiatnre  profiles  ol  the  model 
neutral  atm.  sphere. 


Table  7  illustiates  the  altitude  dependence  of  the  quantities  for 
a  simulation  of  an  H.  release  into  the  atmosphere-  gisen  in 
Figure  2  Thermal  diffusion  becomes  important  on  Is  for  re¬ 
leases  beloss  200-hm  altitude. 

The  chemical  reactions  governing  the  releases  are  taken 
from  /inn  anJ  SuthcrUmJ  [1975],  Chemical  reactions  betsscen 
carbon  dioside  and  the  upper  atmosphere  constituents  are 
negligible.  A  molecular  hsdrogen  release,  bosses  er.  produces 
OH.  O...  and  H.O  hs  the  folloamg  set  of  reactions: 


H.  '  O  •  OH  -  II 


r  ern’/s 


OH  *  O  •  O..  -v  H 
OH  -  II.  •  II .O  i  11 
Oil  -  OH  •  11.0  s  O 


k,  =  5  x  10  "  cm5  s 

A,  -  C6  X  10  *-w  7  cni’/s 

k.;  :  m  10  T  env'/s 
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3G? 
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t  V  Is 


M,  CONCENTRATION 
IOC. .(Number  Density. ern  *) 


(a) 


time  ir»4  i:\y.c 

M  -t.y  0  <2 


H,  CONCENTRATION 
LOG. .(Number  Density, cm  ’) 


(b) 


TIME  -  1J>4  I'SEC 
123  180  -1:3  -62  O  02  123  18G 
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H,  CONCENTRATION 
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[or  (.’()..  Inc  production  and  loss  teims  in  (?|  are  zero.  The 
I li r’. ctt. c'.-*  /’,  and  /..  lor  the  hydrogen  release  are 

/’,  --  (I  /.,  -  A0/r,;.\, 

k»r..r.„  ~  (*,«..  4  *,«,  +  kan,)n. 

" here  !he  subscripts  I  and  2  refer  to  the  H..  and  OH.  respec- 
1 1 s els .  and  rr„  is  the  atomic  oxygen  eoneentration. 

I  he  neutral  gas  expansions  arc  illustrated  as  isopycnic  con¬ 
tour  plots  or  in  terms  of  three-dimensional  surface  plots.  'I  he 
gas  cloud  is  outlined  by  a  surface  on  which  the  concentration  is 
10"!  of  that  of  the  central  cloud 

I  lie  conditions  for  loi  r  simulations  produced  by  the  model 
are  listed  in  Table  3.  Simulation  I  modeled  the  point  release  of 
3  >  10”  molecules  of  CO,  and  H,  at  300-Ism  altitude.  The 
neutral  wind  velocity  was  set  to  zero  Figure  3  illustrates  the 
diffusive  expansion  31.9  s  after  release.  The  CO,  cloud  is  more 
concentrated  but  smaller  in  si/e  than  the  faster  diffusing  H. 
cloud.  The  maximum  OH  concentration  is  0.2  -  of  the  maxi¬ 
mum  II,  concentration.  At  the  early  times  of  expansion  the 
clouds  arc  spherical.  I  igure  4  illustrates  the  II,  and  CO,  gases 
almost  4  min  after  release.  The  increase  in  the  neutral  diffusion 
coefficient  with  altitude  causes  nonspherieal  expansions  at  this 
later  time. 


I  igure  5  illustrates  the  effect  of  nonspherieal  expansion  on 
the  maximum  concentration  inside  the  expanding  clouds.  The 
maximum  concentration  for  spherical  expansions  is  propor¬ 
tional  to  ;  3  •’  as  shown  by  the  dashed  lines  in  the  figure.  1  he 
actual  expansions  deviate  from  the  dashed  lines  for  times 
greater  than  1  min. 

The  buildup  and  decay  of  the  OH  radical  are  illustrated  in 
Figure  6.  I  he  OH  concentration  reaches  a  maximum  value  of 
I  05  X  10' cm  3  at  II  s  after  the  H,  release.  The  OH  concentra¬ 
tion  never  becomes  more  than  0.3'!  of  the  H,  concentration 
for  releases  at  300- Is m  altitude. 

The  results  of  simulation  1  are  compared  to  an  approximate 
analytic  solution  for  diffusion  in  an  exponential  atmosphere. 
The  purpose  of  this  comparison  is  to  provide  an  estimate  for 
conditions  when  the  approximate,  analytic  expression  can  be 
used. 

The  approximate  expression  for  the  gas  concentration  re¬ 
sulting  from  a  point  release  is  given  by 
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I  ig  U  I  xpansion  of  a  vertical  H,  line  source  centered  at  2'0  k in  altitude 
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w  here  \„  is  the  number  of  molecules  released,  />„  and  II,  are 
the  diffusive  coefficient  and  scale  height,  respective!),  of  the 
atmosphere  at  the  point  of  release,  ami  //,  -  17  »ij  is  the 
released  gas  scale  height. 

This  formula  was  derived  bv  Bernhardt  [  1976)  and  has  been 
Used  bv  Bernhardt  cl  at.  [  1  7 5 ] ,  Bernhardt  and  da  Rosa  [1977], 
Mer.ddio  and  I'orhes  [I97S],  and  Anderson  and  Bernhardt 
[  197s]  to  describe  the  diffusive  flow  of  point  releases. 

The  comparison  is  made  for  the  H,  release  at  300-km  alti¬ 
tude  because  the  elVects  of  thermal  diffusion  and  chemical 
reaction  are  minimal  The  two  solutions  are  identical  for  the 
first  40  s  All  54.5  s  after  the  release  the  peak  concentration  of 
the  an.dv  tic  solution  is  larger  than  that  of  the  numerical  solu¬ 
tion  Also,  the  anal) tic  solution  extends  upward  slight!)  more 
than  the  numerical  solution  (f  igure  7). 

\t  later  times  the  difference  between  the  analytic  and  nu- 
meric.d  solutions  becomes  more  pronounced  Most  of  this 
deference  is  in  the  cloud  topside,  figure  8  illustrates  the  nu¬ 
meric, d  and  analytic  solutions  507. X  s  after  release  W  e  find  the 
anal) tic  solution  to  he  a  useful  approximation  to  the  numeri¬ 
cal  solution  for  the  first  10  min  for  a  release  near  300-km 
altitude 
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I  he  lesults  of  simulation  2  are  shown  111  I  igure  9  An  II, 
souKe  moves  veiticallv  with  a  velocity  of  X  r. m  s,  rele.isim,' 
vapor  at  .1  constant  rate  of  6  /  I  O'4  molecules  s  from  230-f  nt 
to  270-kin  altitude.  Since  the  release  starts  at  f.r.s  altitudes  and 
moves  upward,  the  gas  cloud  bulges  at  the  bottom  during  the 
earlv  times  ol  expansion  (figure  9 a),  for  expansion  limes 
greater  than  20  s  the  iinnuniforinilv  of  the  background  atrno- 
plmre  causes  a  bulge  in  the  upper  portion  of  the  cloud.  Mans 
minutes  into  the  expansion  the  gas  loses  all  memory  of  its 
initial  structure  (I  iguie  9c) 

1  he  elicits  ol  neutral  winds  are  demonstrated  in  the  next 
two  simulations  The  neutral  wind  profiles  Used  in  the  simula¬ 
tion  (illustrated  in  I  iguie  10)  are  taker,  from  measurements  at 
Sondre  Stromfjord  Station  ((>7  02'N,  50. 17°VV)  011  December 
1 7,  l*D4.  as  reported  bv  t'ererra  cl  cl  [  197k].  The  strong  shear 
in  tins  w  md  has  a  striking  ellect  on  the  neutral  gas  expansion 

I  igute  I  I  illustrates  the  effect  of  a  nonuniform  neutral  xsind 
on  this  diffusion  of  3  x  I0;:  molecules  of  CO-  released  at  a 
point  l no  km  in  altitude  T  he  neutral  wind  shear  elongates  the 
initially  spherical  gas  expansion. 

The  last  simulation  illustrates  the  effect  of  the  strong  neutral 
wind  shear  on  an  elongated  source.  A  CO,  source  moves  at  453 
to  the  vertical  vv  ith  a  velocity  of  1 .4  kni/s.  releasing  vapor  from 
145-km  to  175-km  altitude  at  a  rate  of  I  X  ICP  molecules  s. 

I  lie  neutral  wind  pattern  changes  the  initially  west-east  gas 
distribution  into  a  north-south  distribution  covering  tens  of 
kilometers  (figure  12). 

7.  St  xtM  xKV  wn  Cosrt.tsioss 

A  iHimerie.il  model  of  neutral  diffusion  has  been  developed 
for  the  study  of  gas  flow  in  a  nonuniform,  chemically  reactive 
environment.  The  three-dimensional  model  permits  realistic 
simulation  of  neutral  gas  expansions  under  a  wide  range  of 
atmospheric  conditions. 

Recently.  Anderson  and  Bernhaidt  [  1 97X]  have  shown  that 
molecular  hydrogen  releases  into  the  equatorial  ionosphere 
will  produce  ionospheric  depletions  which  can  trigger  a  Ray¬ 
leigh- fav  lor  type  gravitational  instability.  They  also  have 
shown  that  H,  releases  will  produce  airglow  at  6300  A  and 
3060  A  due  to  excited  O  and  OH.  Their  study  was,  however, 
limited  to  point  releases  of  H,  using  the  analytic  formulation 
for  gas  expansion  by  Bernhardt  [1976]  Similar  limitations 
were  experienced  by  Bernhardt  cl  al  [  I975J  and  Bernhardt  and 
da  Rosa  [1977]  in  their  studies  of  mid-latitude  ionospheric 
modification.  Research  is  currently  underway  using  the  three- 
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dimensional  neutral  expansion  model  described  here  along 
with  the  existing  ionospheric  models. 

Optical  observations  of  neutral  chemicals  released  into  the 
upper  atmosphere  have  provided  infi'rmution  about  the  wind 
velocities  and  ditVusiou  coefficients  in  this  region.  For  example, 
trimethvl  aluminum  reacts  with  atomic  oxygen  in  the  atmo¬ 
sphere,  leaving  a  chemiluminescent  trail  [Rosenberg  et  a!.. 
1963).  Simulation  of  chemiluminescent  reactions  is  analogous 
to  the  simulation  of  the  molecular  hydrogen  and  atomic  oxv¬ 
gen  reaction  described  here.  Interpretation  of  chemical  tracer 
experiments  may  be  aided  By  u*e  of  the  three-dimensional 
diffusive  expansion  model. 
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C.  GAS  DYNAMICS  OF  VAPOR  RELEASES 


We  have  investigated  the  flow  of  vapors  released  into  the  upper  atmo¬ 
sphere.  Our  previous  work  (see  Section  B)  presented  a  three-dimensional 
model  of  neutral  gas  diffusion.  This  work  considered  the  transition  from 
a  collisionless  expansion  (immediately  after  release)  to  a  collision- 
dominated  (i.e.,  diffusive)  flow.  This  study  was  necessary  to  determine 

the  stopping  distance  of  vapors  released  into  the  tenuous  upper  atmosphere 

5 

where  the  mean  free  path  is  approximately  10  meters  at  400  km  altitude. 

The  details  of  our  new  model  are  presented  in  the  following  paper 
which  was  published  in  the  Journal  of  Geophysical  Research.  Briefly,  the 
vapor  flow  process  is  described  by  the  Boltzmann  equation  with  a  collision 
term  rep. jsenting  relaxation  of  the  injected  vapor  to  the  kinetic  state  of 
the  background  atmosphere.  Simulations  show  significant  differences 
between  supersonic  and  subsonic  vapor  releases.  Lower  speed  releases  take 
a  longer  time  to  become  col  1 isionally  dominated.  High  speed  flows  expe¬ 
rience  large  heating  due  to  collisions  with  the  background  atmosphere. 

The  new  model  permits  prediction  of  the  exapnsion  of  vapors  that  are 
released  into  the  upper  atmosphere  from  supersonic  vehicles. 
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High-Altitude  Gas  Releases:  Transition  From  Collisionless 
Flow  to  Diffusive  Flow  in  a  Nonunil'orm  Atmosphere 

Pai  i  A.  Hi  kmi  xrdi 

Ru.hu..ien.  t'  I  uhui.iluri .  StunJ.-rd  l  'III  i  f  -III.  SlunfurJ.  Califurmu  WT'A 

Vapors,  released  at  high  altitudes  from  rockets  or  the  Space  Shuttle,  m;i)  travel  mans  hundreds  ol 
kilometers  before  coming:  to  rest  An  understanding  id  the  elfecls  produced  In  these  releases  requites 
knowledge  of  the  vapor  flow.  In  the  tenuous  upper  atrnosphere,  the  injected  gases  c V P'.  d  Irvun  a 
colhsionless  to  a  collision-dominated  state.  In  this  paper,  this  process  is  described  bv  the  If  jlt/maun 
equation  w  ilh  the  krviok  collision  integral  I  lie  tlieor)  is  applied  to  supersonic  releases  intsi  a  nonunifom 
aliuosphere.  The  mode!  predicts  elongation  and  heating  of  the  vapor  trails  due  to  Collisions  I  he  theory 
may  be  applied  to  a  wide  range  of  experiments  whose  objectives  include  the  formation  of  ionosphere 
depletions,  the  triggering  of  instabilities,  the  generation  of  gravity  waves,  and  the  injcciivm  of  neutrals  as 
luminous  tracers. 


I.  Introduction 

The  flow  of  vapors  released  into  the  upper  atmosphere  is 
governed  by  three  types  of  expansions  During  the  early  stages, 
collisions  between  the  expanding  vapor  molecules  dominate 
the  flow.  This  expansion  is  called  self-continuum  How.  As  the 
expanding  cloud  becomes  tenuous,  self-collisions  become  rare 
and  in  a  vaccutrt  the  expansion  is  free  streaming  or  collision- 
less.  At  some  time  in  the  expansion,  collisions  with  the  back¬ 
ground  atmosphere  become  important.  When  interaction  with 
the  background  atmosphere  controls  the  flow,  a  state  of  dill  it- 
sive  expansion  exists. 

The  nature  of  the  expansion  depends  on  (he  amount  of 
material  released  and  on  the  par’icle  concentration  of  the 
background  atmosphere.  In  a  der.,.e  atmosphere,  the  expan¬ 
sion  interacts  with  the  background  atmosphere  to  form  R.tnk- 
ine-llugoniot  shock  waves.  These  shock  waves  separate  the 
background  atmosphere  from  the  core  of  the  expansion.  When 
releases  are  made  at  high  altitudes,  i.e  ,  at  reduced  atmospheric 
concentrations,  the  shockwaves  are  broader,  and  the  transition 
between  expansion  core  and  the  background  becomes  levs  well 
defined  Kurt  her  reduction  of  the  background  concentration 
eliminates  any  shock  waves.  K.xperimental  verification  of  tbeve 
flow  regimes  is  presented  by  Muni;  el  a I  [I97U[. 

All  gas  expansions  initially  are  governed  by  continuum  flow , 
Such  How  in  a  vaccunt  has  been  described  by  Freeman  urn! 
(irumiy  [I9oS]  and  Freeman  ami  Thomas  (I9(>9).  Continuum 
How  in  a  raritied  atmosphere  lias  been  described  by  llohuiy 
(|9b9j.  Trade  and  Tristram  [1972],  and  llojjnum  and  llesi 
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Ihe  transition  from  self-continuum  How  to  dtlltisivc  How 
has  been  examined  by  liman  and  Baxter  [I977j.  In  their 
treatment,  they  assume  that  the  gas  velocity  distribution  does 
no!  deviate  signtlicanlly  from  a  Maxwellian  distribution.  I  he 
How  is  collision  dominated,  either  by  self-collisions  or  by 
viuh'tons  with  the  background  atmosphere,  throughout  the 
expansion  I  his  tv  pe  of  solution  is  valid  in  atmospheres  denser 
than  those  considered  lit  this  paper. 

In  a  low -density  background,  the  transition  front  a  self- 
v 1 1 ;  tin  it  uni  e*  pan  sum  to  a  dill  u  six  e  expansion  max  be  divided 
into  two  expansion  sequences  Ihe  first  sequence  consists  of 
transition  liom  sell-conlinttiini  to  colhsionless  How  [Draper 
and  Dili.  I  Wtj  I  he  v.Nond  sequence  involves  the  transition  of 
kollisionlcss  How  to  dilhisixe  How.  I  his  latter  transition  has 
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been  modeled  by  Hrauk  and  Hamel  [1972],  and  Baum  (19 '2, 
197-1]  for  interaction  with  a  uniform,  chemically  nonre.tcitve 
atmosphere.  At  high  altitudes  the  self-eonliriuum  flow  is  un¬ 
coupled  from  the  diffusive  flow  by  the  intervening collixionle's 
expansion. 

Th  e  work  presented  here  extends  the  efforts  of  Baum  [197  s. 
197-1]  to  include  the  effects  ol  a  nonuniform  background  atmo¬ 
sphere  and  a  velocity  dependent  collision  frequency.  The  un¬ 
steady  expansion  of  gas  releases  is  analyzed  using  gas  kinetic 
theory  The  Boltzmann  equation  with  the  Krook  collision 
integial  is  solved  numerically  At  late  times  (after  many  colli¬ 
sions  with  the  background  atmosphere)  the  solution  is  identi¬ 
cal  with  the  one  given  by  diffusion  in  a  notuinilorm  atmo¬ 
sphere  [Bernhardt.  I979| 

This  analysis  is  motivated  by  a  need  for  a  description  of 
high-altitude  releases  from  rockets  or  from  the  Space  Shuttle 
Bernhardt  [I97(>],  Mer.ddlo  and  lathe-.  [I9,'!vj,  and  Andersm 
amt  Bernhardt  [I97S]  have  show  a  that  the  release  into  the  f 
layer  ionosphere  of  certain  iieulial  gases  (such  as  CO;.  Il-O. 
and  II;)  will  produce  plasma  depletions,  plasma  temperature 
enhancements,  and  airglow  emissions  and  may  trigger  in¬ 
stabilities.  Proper  interpretation  of  the  results  from  such  ex¬ 
periments  requires  modeling  of  the  How  of  the  released  gas 
This  is  especially  true  for  the  supersonic  releases  made  from 
the  Space  Shuttle. 

This  research  was  also  motivated  by  the  observations  of 
Sehunk  [I97S]  that  a  simple  diffusive  model  may  not  ade¬ 
quately  describe  the  floxv  resulting  from  gas  releases  Specifi¬ 
cally.  Sehunk  states  that  the  effects  of  nonlinear  acceleration, 
viscous  stress,  and  velocity -dependent  collision  frequencies 
may  significantly  retard  the  expansion  of  the  gases  at  high 
altitudes.  The  solution  to  the  Boltzman's  equation  incorpo¬ 
rates  all  of  these  effects. 

2.  Avm  vnr  Di  sr  riim  ion 

In  kinetic  theory,  the  stale  of  a  gas  is  completely  described 
by  a  distribution  function  /(.x.  v,  I),  where  x  is  the  position 
vector,  v  is  the  velocity  vector,  and  t  is  time.  The  first  three 
velocity  moments  of  this  function  give  the  common  physical 
properties  of  the  gas 

/r/'v  -  n(x.  I) 

J  \JtPv  -  u(x.  r)u(x,  r) 
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where  n  is  eoneentr.UH'it,  u  is  hulk  velocity ,  anil  7 ( x,  /)  is 
temperature 

I  he  evolution  of  the  gas  distribution  function  is  described 
bv  the  Boltzmann  equation 

-L  *  v-VJ  +  F-VJ  -  (ef/Pt),  (1) 


expanded  to  the  degree  that  the  self-collision  frequency  tv, )  is 
essentially  zero  1  he  atmosphere  prov  ides  a  large  momentum 
and  energy  sink  for  the  released  molecules  [Baum.  1973].  Also, 
the  high-velocity  molecules  are  primarily  moving  into  undis¬ 
turbed  regions  of  the  atmosphere.  Consequently,  the  back¬ 
ground  atmosphere  is  assumed  to  be  in  an  unperturbed  state 
throughout  the  expansion.  Under  these  assumptions,  (4)  is 
w  rrtten  as 

4^  ■+  v-VJ  -t  F-VJ  -  •'„(*,  v)(nr/>0  -  /) 


where  V,  and  V,  are  the  gradient  operators  in  x  and  v  space, 
respectively;  F  is  an  external  acceleration  vector,  and  (<  f/cl), 
is  the  change  in  the  distribution  due  to  collisions.  In  this  note, 
the  krook  collisional  term  is  used, 

-  »■(«•>  -  /)  (2) 

where  r(x,  v)  is  the  collision  frequency.  (3  is  the  Maxwell 
distribution  given  by 

(  in  V 1 

*(V)=  CXP 


+  &(»*Me)*[x  -  S(t)l  (6) 

where  the  subscript  e  has  been  dropped  for  convenience. 
Coordinates  x.  y ,  and  :  denote  the  cast,  north,  and  vertical 
directions,  respectively.  The  effects  of  gravity  can  be  repre¬ 
sented  by  F  =  -g0 a,,  where  £„  is  the  gravitational  acceleration, 
and  a,  is  the  unit  vector  along  the  vertical. 

The  collision  frequency  must  be  chosen  to  represent  mo¬ 
mentum-transfer  collisions  between  an  exhaust  particle  with 
velocity  v  and  the  background  atmosphere.  The  number  of 
collisions  per  unit  time  is  written  as 
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in  is  the  molecular  mass,  and  k  is  the  Bolt/mann  constant.  The 
Krook  collision  term  represents  relaxation  of  the  distribution 
function  to  a  Maxwellian  distribution  characterized  by  a  tem¬ 
perature  T  and  a  velocity  u  [Bhalnagar  el  at..  1954], 
Simulation  of  gas  releases  involves  two  distribution  func¬ 
tions:  /,  the  released  gas  (or  exhaust)  function  and  /„  the 
atmospheric  distribution  function.  The  coupled  Boltzmann 
equations  describing  the  interaction  between  these  functions 
are 


-~y  +  v- V,/,  +  F •  VJ,  -  -  /,)  +  -  !<) 

+  ~  Su)]  (4) 

+  V-Vj0  t  F-V,./a  =  l'a(/ljOa  -  /„)  +  I'ofl'hO,  -  fa) 

r'  (f) 


wlrcrc  the  subscripts  e  and  a  refer  to  the  released  gas  and  to  the 
atmosphere,  respectively;  and  i'0,  represent  collisions  be¬ 
tween  the  exhaust  and  the  atmosphere;  (2,(0  is  the  source  flow 
rale.  <>.( v )  is  the  initial  source  distribution  function;  and  J[x  — 
S(r)j  represents  a  moving  point  source  with  position  vector 
S(r)  -  Sxa,  3-  ,S\ov  +  ,V,a,.  The  values  used  for  u  and  T  in  the 
lunction  are  indicated  by  the  subscripts  e,  a,  and  s  on  the 
distribution  function  d>.  The  molecular  masses  used  in  (4)  and 
(5)  are  m,  arid  ma,  respectively.  The  coupling  between  these 
equations  occurs  through  the  collision  frequencies  (»■,„  and 
»•„,).  which  are  concentration,  velocity,  and  temperature  de¬ 
pendent.  Production  and  loss  terms  can  be  added  to  (4)  and  (5) 
tl  chennc.il  reactions  are  important. 

All  of  the  gas  flows  described  in  the  introduction  can  be 
modeled  by  solution  of  (4)  and  (5).  For  example,  if  all  the 
collision  frequencies  arc  set  to  zero,  collisionless  flow  occurs.  If 
e,  »  ,  taking  the  velocity  moments  of  (4)  yields  the  equa¬ 

tion',  governing  self-continuum  flow.  Similarly,  neglecting  r, 
and  finding  velocity  moments  can  provide  the  equations  de¬ 
scribing  dillusive  flow. 

In  this  paper,  the  discussion  is  limited  to  the  transition  from 
collisionlcss  to  diffusive  flow  The  gas  is  assumed  to  have 


v0(x,  v)  ,  zi0(x)<j,0tVo(v ) 

nif  t  nia 

[Banks  and Koekans.  1973]  where  m,  and  ma  arc  the  masses  of 
the  exhaust  and  atmospheric  particles,  respectively,  n„(x)  is  the 
atmospheric  concentration.  ir,a  is  the  collision  cross  section, 
and  tf0(v)  is  the  average,  relative  speed  between  the  atmo¬ 
spheric  molecules  and  the  exhaust  molecule  with  velocity  v. 

The  atmospheric  concentration  is  assumed  to  vary  ex¬ 
ponentially  with  altitude; 

rto(c)  =  «oo  exp  [(z„  -  r )///«] 

where  iiaj  is  the  concentration  at  altitude  z0,  and  Ha  =  kT0/ 
ntjg  is  the  scale  height.  The  values  for  temperature  T0  and 
average  molecular  mass  w„  are  found  from  the  MS1S  neutral- 
atmosphere  model  described  by  Hedin  el  at.  [1977o,  6],  The 
altitude  dependence  should  be  included  if  the  mean-free  path  is 
greater  than  the  scale  height  of  the  atmosphere. 

The  collision  cross  section  is  =  ir (r,  +  r„y.  The  exhaust 
and  atmosphere  molecules  are  taken  to  be  rigid  spheres  with 
radii  r,  and  ra.  respectively.  With  this  model,  the  collision  cross 
section  is  independent  of  velocity.  The  molecular  radii  for  air, 
COj,  and  H,  are  0.214,  0.20,  and  0.145  nm,  respectively  [Uir- 
schf elder  el  al..  1954], 

7  he  differential  speed  between  the  injected  gas  molecules 
and  the  atmosphere  is  found  by  averaging  the  relative  speeds 
over  the  atmospheric  velocity  distribution 

uffv)  =  J  (»')</ V 

w here g  -  ]v  -  w|  is  the  relative  speed  between  an  exhaust 
molecule  with  velocity  v  and  an  atmospheric  molecule  with 
velocity  w.  Substitution  of  the  Maxwell  distribution  given  by 
(3)  and  integration  yield 
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Fig  I.  lncrc;i>c  in  the  molecular  collision  Deque  re  >  wilh  speed 
The  exact  expression  involves  exponentials  and  error  fur.*itonv  The 
approximate  formula,  which  only  requires  the  computation  of  a 
square  root,  is  used  in  the  modeling. 


where  ca  =  {2k  Ta / nia )‘ 2  is  the  atmospheric  speed  of  sound  and 
i  v  -  Uu  |  is  the  relative  speed  between  the  exhaust  particle  and 
the  atmosphere  [\f uckenfuss.  1 962 J.  l  or  |  v  -  ,  <\  «  I,  the 

relative  speed  becomes 

2 

t»a  ~  j- 2  Ca  -  (&A  I  a/~  tHa  ) 

7T 

Thus  al  loss  speed  the  average  relative  speed  (and  also  the 
collision  frequency)  is  independent  of  exhaust  velocity.  If  j  v  - 
u j ■'  /<•„  »  I ,  the  average  relatis e  speed  becomes  i ,  ,(v )  =•  jv 
u„| .  At  high  speeds  the  thermal  motion  of  the  atmosphere  has 
little  effect  on  the  average  relative  velocity. 

liquation  (7)  may  be  approximated  by  a  hyperbola  to  sim¬ 
plify  numerical  calculation 


The  maximum  error  caused  by  this  hyperbolic  approximation 
is  less  than  3r< .  The  expression  for  collision  frequency  becomes 


l  L  1 

ft1  1. i4  i.V 
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I  !.*  '  Aluiiulc  .mil  Mil.tr  cxclc  vlcpcrdcncc  of  the  mc.tr.  I;cc  p.ith  for 
<  ( >,  hi  ihi-  .timosphcic 


where  i„  -  (2ia«i„n,1()<j,„  )/[ir‘  ’{in,  +  m„)].  The  exact  and  ap¬ 
proximate  computations  for  the  velocity  dependence  of  the 
collision  frequency  are  compared  in  f  igure  1.  1  he  rapid  in¬ 
crease  in  the  collision  frequency  for  velocities  greater  than  the 
speed  of  sound  is  illustrated. 

Another  parameter  of  interest  is  the  mean  free  path  of  each 
exhaust  particle  as  defined  by 

A(x,  v)  -  |  v  -  u„ | /l'„(X,  v) 
for  low -speed  particles, 

A, *.  v)  , 

"'s'lofllK.CI. 

which  vanishes  at  |v  -  u^l  -  0.  Tor  high-speed  (supersonic) 
particles,  the  mean  free  path  lakes  the  limiting  value 


A(x.  v)  • 


m,  -f  m„ 
m0na(%  )ara 


Representative  profiles  of  the  high-speed  values  of  mean  free 
path  for  CO,  are  shown  as  a  function  of  solar  cycles  in  figure 
2.  These  profiles  arc  based  on  the  MS1S  model  atmosphere.  A 
strong  altitude  and  solar  cycle  dependence  is  evident. 

1  he  formulation  represented  by  (6)  is  only  valid  in  low- 
density  regions  where  the  released  gas  w  ill  attain  acollisionless 
state  before  the  flow  is  cominated  by  diffusion.  The  regions, 
which  permit  collisionless  flow,  are  found  by  comparing  the 
equal-mass  radius  of  the  released  gas  to  the  atmospheric  mean 
free  path  (e  g.,  Baxter  and  Linton.  1977;  Linson  and  Baxter. 
1977],  The  equal-mass  radius  is  the  characteristic  length  of  the 
cloud  when  it  has  expanded  to  a  volume  which  contains  an 
amount  of  atmospheric  mass  equal  to  the  mass  of  the  released 
gas.  This  radius  is  given  by 


where  is  the  number  of  exhaust  molecules  released  and  rt „  is 
the  atmospheric  concentration. 

If  this  length  is  less  than  3  mean  free  paths,  the  released  gas 
will  expand  in  a  collisionless  state  for  at  least  1.5  mean  colli¬ 
sion  times  before  significant  interaction  with  the  background 
atmosphere  occurs.  Consequently,  (6)  can  be  applied  if «,/ A  < 
3.  This  criterion  is  illustrated  in  figure  3  for  a  CO,  injection 
into  the  MSIS  model  atmosphere  during  equinox,  1977. 

The  procedures  described  in  this  paper  may  be  applied  for 
conditions  of  release  above  and  to  the  left  of  the  curve  in 
figure  3.  Releases  corresponding  to  conditions  below  the 


l'i>:.  3  Conditions  when-  different  expansion  concepts  apply  The 
theory  presented  <n  the  text  is  applicable  to  releases  for  conditions  to 
llic  Id t  of  ami  abase  the  solid  curse. 
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cunc  should  use  the  continuum  flow  equations  such  as  de¬ 
scribed  by  Hauer  and  Liman  [  1977).  1  he  validity  of  the  criteria 
described  above  is  only  approximate.  Changes  in  the  initial 
velocity  and  shape  of  the  injected  gas  affect  these  criteria. 

3.  M  l  lllOK  Ut  Sot  l.TION 

I  quatioil  (6)  is  a  partial  differential  equation  in  seven  vari¬ 
ables  r,  r,  t\,  t\,  t,,  and  i.  An  analytical  solution  for  a 
stationary  source  [S(r)  -  0]  in  a  uniform  atmosphere  (i-Q  - 
const)  has  been  found  by  Bantu  [1973,  1974]  This  analytic- 
solution  involves  a  three-dimensional  convolution  in  space  as 
well  as  an  integration  over  time.  Baum's  solution  is  not  readily 
extendable  to  the  analysis  of  releases  from  a  moving  source  in 
a  nonuniform  atmosphere  with  a  velocity-dependent  collision 
frequency. 

If  the  flow  were  nearly  Maxwellian,  velocity  moment  equa¬ 
tions  could  be  used  to  describe  the  solution.  Collisionless  flow  , 
however,  is  highly  non-Maxvvcllian,  and  truncation  to  a  finite 
set  of  velocity  moment  equations  would  introduce  large  errors 
in  the  results.  For  this  reason,  a  numerical  solution  to  (6)  is 
formulated. 

By  transformation  of  variables  and  integration  over  velocity 
space.  (6)  is  converted  into  an  integral  equation  suitable  for  a 
recursive,  numerical  solution.  Using  the  transformations 


w  here 

G;(x.  v,  r,)  =  exp 


'"(*  -  v'.  +  t  -  i.)  C'u(v  -  Fr,)  dt,  (11) 
|j-  J  '  Mx  -  vr,  +  F/jV2,  v  -  Fr2)  r/tjj 


/(a,  c,  r )  =  0  r<0 


The  distribution  function  given  by  ( 1 1 )  is  multiplied  by  1 ,  v. 
and  j  v j 1  and  integrated  over  velocity  space  to  yield  concentra¬ 
tion,  bulk  velocity,  and  temperature,  respectively.  The  result  is 
written  compactly  as 


n(x.  f) 


I 

u(x. /) 

|  u| 2  +  }kT(x,  »)/m,  . 


I. 


C3[x.  S(r-  r,),/l)-^('-J-l)-d,,[U(r,)) 
*1 


/  j  G3(x,  x',  V((j) 

J0  J  -  a> 


r  1  1 

U(/,)  j 

JU(/,)i2  j 


dh 


x'  =  x  -  vr  + 


Fr2 


and 


v’  -  v  -  Fr 


«(x’.r  -  r,) 


d>«[V(r, )] 


r  1 

i  V(r, ) 

L  I V(r,)|J 


t/’x'  dt,  (12) 


F.quation  (6)  becomes 

-  eu(p.  q)[u(p,  t)'3j(q)  -  /]  +  (MOvUpWP  ~  S(r)]  (9) 

vv  here 

p(r)  =  x  +  v'r  +  Fr2/2 
q(r)  =  v'  +  Fr 
\\  ith  the  substitution 

f  =  g  exp  j-  f'VolP(ti).q(t,)]d/ 

Fquation  (9)  becomes 


-f-  =  G,(x’,  v',  r )<^,( / )o[p  -  S(r)) 

<  r 

+  G,(x’,  v’.  r)i'„(p,  q)u(p.  lk'»„(q)  (10) 


w  here 


vhere 


U(r,)  =  [x*  S(r  -  r,))/r,  -  Fr,/2 
V(r ,)  =  (x  -  x’)/r,  -  Fr,/2 
x’  =  x  -  vr,  +  Fr,2/2 
Gjfx.S.r,)  =  Gj(x. (x  -  S)/r,  +  Fr,/2.r,) 

The  first  equation  in  (12)  is  a  Fredholm-Volterra  equation 
of  the  second  kind.  There  exist  a  variety  of  methods  for  its 
solution  [e.g..  Baker,  1977],  Quadrature  techniques  are  chosen 
for  simplicity.  The  first  term  on  the  right  side  of  (12)  is  eval¬ 
uated  using  Gaussian  quadrature.  This  is  the  dominant  term 
for  time  less  than  I/p.  Early  lime  solutions  arc  introduced  into 
the  second  term  on  the  right  side  of  (12).  The  solution  then 
proceeds  recursively.  After  the  concentration  n(x,  r)  is  found, 
the  velocity  and  temperature  are  found  directly  from  the  sec¬ 
ond  and  third  equations  in  (12). 

In  certain  cases,  the  integral  representation  of  G ,  may  be 
computed  analytically.  If  the  external  acceleration  term  (i.c., 
gravity)  is  dropped,  G,  becomes 


v'.  and  g  yield 


'(x.  v.  I) 


G3(x,  S,  i,)  =  exp 

w  here 

j  ( c  x  -  s\  1 

<uIpf/, ).  q(/,)l  dt,j 

and  making  substitutions  for  x’. 

/  AT )  -- 

exp  [(.-  -  -V,  V/fd)  ~  1 
(c  -  AT)  //a 

If  both  external  forces 

and  atmospheric  gradients 

ii  )'.'>.(  v  -  Fl,) 

glccted,  /•(:,  AT)  -  1. 

{x  vt,*  -  S(r-  /,  )J  </r, 
t  J  G,(x.  v.  l,)  r,  (x  -  vi,  4  — v  -  FT,) 


In  the  following  discussions.  F  is  set  to  zero.  This  is  justified 
because  gravity  is  negligible  for  the  short  limes  considered  in 
the  examples.  Also,  releases  at  orbital  velocity  do  not  see  the 
ellects  of  gravity  in  a  flat  earth  geometry.  A  more  detailed 
analysis  may  have  to  include  the  effects  of  a  curved  (versus 
flat)  earth  as  well  as  gravity. 
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Fig.  4.  Simulation  of  a  low-tcmperaturc  ( 300  K )  spherical  expansion  of  COj.  The  concentration  and  temperature  of  the 
background  pas  are  1  X  10"  m  3  and  870  K,  respective!)  The  loss-speed  atmospheric  collision  frequency  (t)  is  70.8  s.  At 
early  times  (/  <  r),  the  expansion  is  nearly  collisionless.  At  later  times,  the  concentration,  velocity .  and  temperature  relax  to 
their  collision-dominated  limits. 


4.  Physicai.  Interpretation 

The  released  gas  concentration  given  by  (12)  can  be  written 
as 

n(x,  I)  =  ;t0(x,  I)  +  n,(x,  I)  (14) 

The  first  term.  njx.  t),  denotes  (he  part  of  the  solution  which  is 
in  a  colli  si  on  1-ess  state.  Written  out  in  full,  this  term  is 


This  function  is  a  convolution  of  the  previous  gas  concentra¬ 
tion  with  a  function  of  the  background  atmosphere.  It  is  the 
only  term  remaining  after  many  collision  times.  The  function 
n;(x,  i)  asymptotically  approaches  the  diffusion  solution  for 
long  times  [Baum,  1973]  and  is  called  the  forced  response  to 
the  gas  injection.  In  this  example,  the  forcing  function  is  the 
particle  distribution  of  the  background  atmosphere. 


n,lx.  i )  = 


[  G'j(x. 

»  o 


S(t  -  r,).  i,) 


(?.('  "  f.) 


0. 


x  -  S (t 


('  -  ».)]  ,r 

7 - V1' 


(15) 


where  F  is  set  to  zero. 

This  is  the  only  remaining  term  for  a  vanishing  collision 
frequency .  The  function  C,  (see  (12)  and  (13))  represents 
dumping  of  the  collisionless  expansion  by  interaction  ssith  the 
background  atmosphere.  The  function  «„(x,  /)  is  excited  only 
h>  the  imti.il-g.is  source  distiibulion  function  and  is  called  the 
natural  response 

The  second  term  in  (14)  is  the  part  of  the  solution  governed 
by  the  background  environment. 
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5.  GasSih  rci  Ri  presentation 

Both  point  sources  and  line  sources  of  released  gases  are  of 
interest.  The  release  of  gas  vapors  from  explosively  detonated 
canisters  is  an  example  of  a  point  source  release.  The  line 
source  can  be  represented  by  exhaust  released  from  a  rocket  or 
from  the  Space  Shuttle.  Only  these  sources  will  be  considered 
here. 

The  initial  velocity  of  molecules  released  from  any  source  is 
the  sum  of  the  source  velocity  vector  plus  the  initial  velocity 
vector  of  the  molecules  relative  to  the  source: 

u,  -  v,  4  w,  (17) 

where  u,  is  the  hulk  velocity  relative  to  a  point  in  space,  v,  is 
the  source  velocity  relative  to  a  point  in  space,  and  w,  is  the 
hulk  velocity  relative  to  the  source. 

A  point  source  is  specified  by  the  function  (?,(/)  * 

"here  Q„  is  the  number  of  molecules  released.  T  he  position 
vector  S0,  velocity  vectors  u„.  v,,  and  w„.  and  temperature  T, 
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J  ig.  5.  Simulation  of  a  high-tempcraturc  (?400  K  )  spherical  expulsion  of  CO,  The  background  gas  is  identical  to  the 
one  for  1  tgurc  4  (r  -  70.8  s).  The  large  initial  temperature  at  the  time  of  release  produces  a  rapid  expansion.  At  the  edge  of 
the  expansion,  the  velocity  is  retarded  and  the  temperature  is  desated  above  the  ambient  temperature. 


complete  the  description  of  the  source.  Substitution  of  these 
quantities  into  (15)  yields 

«o(x.  O  -  C3(x.  S„,  l)  &  o.  (IS) 

where  C,  ami  d>>  have  been  previously  defined  in  (?)  and  (12). 
respectively . 

'I  he  line  source  is  assumed  to  be  moving  at  a  constant 
velocity  and  to  be  releasing  molecules  at  a  constant  rate.  The 
position  of  the  source  is  given  by 

S(r)  -  v,r  +  S„ 

where  So  is  the  location  of  the  start  of  the  release  at  time  r  =  0 
1  he  molecules  are  assumed  to  be  released  at  a  rate  of 
molecules  per  second.  The  duration  of  the  release  is  r0  s.  The 
velocity  and  temperature  of  the  molecules  leaving  the  source- 


are  w,  and  T„  respectively.  Using  (IS),  the  natural  term  in  the 
solution  of  the  Bolt/mann  equation  is 

rto(X.  r )  -s  j  C,[X.  v,(r  —  /()  T  So.  b)~  T  o* 

l,nm  o.i  -  i0  *i 

•  +  v.d  -///,)  L,  (19) 

L  «i  J 

6.  Spiiikkm  Hxpanxion 

It  is  instructive  to  consider  the  flow  resulting  from  the  point 
release  of  vapor  into  a  uniform  atmosphere.  The  results  for  the 
spherical  expansion  can  be  compared  with  analytic  solutions 
for  collisionless  and  diffusive  expansions.  The  concentration, 
radial  velocity .  and  temperature  for  the  collisionkss  expansion 
from  a  point  release  are  given  by 


TABI.I:  I.  Simulated  Vapor  I  lows 
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100  0  km  CUBE  CENTERED  AT  400.0  km  ALTITUDE 
TIME  AFTER  RELEASE  2.5  sec 
CO,  CONCENTRATION  AT  SURFACE:  9.2[8]  cm'3 


CO,  CONCENTRATION  AT  SURFACE:  2.1(7]  cm’* 


160.0  km  CUBE  CENTERED  AT  400.0  km  ALTITUDE 
TIME  AFTER  RELEASE:  12.3  sec 

CO,  CONCENTRATION  AT  SURFACE:  9.3(0)  cm*5 


CO,  CONCENTRATION 
LOG, .(Number  Density  .cm-*) 


(d) 


I  ig  7.  Formation  of  ;i  cometlike  tail  on  a  high- velocity  point  release,  moving  hori/omally.  Scattering  of  the  injected 
molecules  cause  the  elongation  of  the  originally  spherical  sapor  cloud 


leased  into  a  uniform  atmosphere  ssjth  a  concentration  of  I  X 
III"  in  \  a  temperature  of  X60  K.  and  an  average  molecular 
mass  of  1X2  amu.  I  he  mean  collision  time  (r)  for  this  atmo¬ 
sphere  is  70S  s.  In  all  figures,  the  radial  dimension  (y> )  is 
divided  hy  the  lime  after  release  (r)  for  display  convenience. 

I  igure  4  presents  the  results  for  a  point  source  of  CO,  with  a 
loss  initial  temperature  of  .'00  K  At  lime  r/2,  the  Ross  is 
ne.irh.  collisionless.  The  concentration,  radial  velocity,  and 
temperature  are  close  to  their  collisionless  limits.  At  time  4r, 
the  s elouty  profile  is  approaching  Us  diffusive  limit  of p/l,  and 
the  g.is  temperature  has  been  elevated  to  700  K  approaching 
the  background  temperature  of  V,0  K 


A  high-temperature  (2400  k)  release  produces  a  much  more 
rapid  expansion  (Figure  5).  At  r/2  (he  gas  is  in  a  transition 
between  collisionless  and  diffusive  flow.  At  4r,  the  core  of  the 
vapor  is  in  a  diffusive  slate.  The  outer  boundary  of  the  gas  ball, 
however,  has  a  radial  velocity  below  the  diffusive  limit  (Figure 
5 b)  and  a  temperature  higher  than  the  background  temper¬ 
ature  (Figure  5c).  These  effects  on  the  surface  of  the  expanding 
sphere  represent  dissipation  of  the  initial  expansion  ci  'ey  of 
the  gas 

These  tw  o  examples  of  gas  expansion  provide  a  ..heck  v  'he 
simulation.  Both  the  low-  and  high-speed  flows  arc  initially 
collisionless,  relaxing  to  nearly  diffusive  flow  after  many  colb- 
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sion  times.  1  he  calculated  Hows  in  both  eases  are  retarded 
ss  hen  compared  with  purely  ditfusive  expansions  Ibis  retarda¬ 
tion  has  been  explained  in  teims  ol  nonlinear  aeeeleration, 
xisexHis  stress,  and  velocity  dependent  collision  frequency  by 
Sjlui'ik  [l‘>7Sj,  all  ol'  which  are  built  into  the  present  model. 

7.  A  niso  i  Kiii'ii  I  xi'xNxttixs 

In  this  section,  fixe  sapor  releases  are  presented,  illustrating 
the  effects  of  a  nonunifortn  atmosphere,  of  a  large  initial 
xeloeitx,  and  of  variations  in  the  released  species.  The  condi- 


100.0  km  CUBE  CENTERED  AT  400.0  km  ALTITUDE 
TIME  AFTER  RELEASE:  10.2  sec 

CO,  CONCENTRATION  AT  SURFACE:  1.4[8]  cm'1 

(a) 


lions  for  the  live  releases  are  listed  in  Table  I.  The  court 
temperature  is  070  K  in  all  eases.  All  of  the  releases  take  p!..c 
at  400-Kin  altitude  in  the  model  atmospheres  given  in  Table  2 
The  atmospheric  wind  is  assumed  to  be  negligible  in  com 
purison  with  the  supersonic  velocities  of  the  expansions 
The  first  example  is  a  point  release  of  CO,  vapor  The  ga 
center  of  mass  is  initially  at  rest.  Because  of  the  nonur.iforr 
atmosphere,  the  expansions  are  different  from  the  ones  pre 
settled  in  the  previous  section.  The  density,  velocity,  and  ten 
perature  of  the  vapor  2r  s  after  release  are  illustrated  m  FTgur 
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100.0  km  CL.' Si:  CENTERED  AT  400.0  km  ALTITUDE 
TIME  AFTER  RELEASE:  10.2  sec 

COa  CONCENTRATION  AT  SCRr’ACE:  4.0(7)  cm  3 

(a) 
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f  <g  C( cxK.i lit!  trail  piodiucd  bv  the  Space  Sin: tile  mov  me  .1;  '  ’  km  s.  fume  C.s  ene.nes  dow  rmurd  with  an  exhaust 
v c!t»v'  1 1 v  ol  }  I  km  x  I  he  engines  hum  for  IU  s  releasing  U>':  molecules  per  second 


h  (At  tin*  point  ol  telcase.  tin*  mean  collision  1 1 1 c  is  7  T~ 
JV.-s  v  )  I  Jlic  expansion  is  sliyhtlv  nonspherical.  ax  indicated  bv 
the  larger  distance  between  contour  lines  above  than  below  the 
release  point  (1  mure  (hi)  At  a  given  distance  from  the  teleasc 
point,  the  expansion  speed  increases  with  altitude  (1  igure  (>/*) 

1  he  vertical  temperature  profile  through  the  center  of  the 
release  1  |  igure  (\  )  shows  greater  eo(lisn>na(  heating  below  the 
release  pmnt  he. ause  of  the  denser  atmosphere  at  low  alti¬ 
tude'  1  xcept  foi  the  small  anisotrupv  introduced  h\  the  non- 
lo’eir  atmosphere,  the  expansion  is  qualitative!)  the  same  as 
st;ovv  n  i n  sev.  lion  <v 

I  lie  second  example  illustrates  the  elongation  ol  a  point 
'"■'M'  when  it  is  released  h-»ri/or,t.ill\  jt  M.^h  I]  into  the 


atmoxphere  during  a  period  of  average  solar  activity.  The 
injection  velocity  is  directed  southward  Successive  time 
Iramcs  ol  the  gas  expansion  are  shown  in  f  igure  7.  Collisions 
with  the  background  gas  molecules  cause  a  comethke  tad  to 
form  on  the  cloud  After  12.3  s  the  vapor  trail  extends  over 
Id0*km  distance  (f  igure  7c/),  and  the  bulk  of  the  vapor  is  still 
in  a  collfsionless  state. 

If  the  vapor  is  injected  with  a  downward  velocity  com¬ 
ponent.  it  will  be  expected  to  be  rapidly  stopped  as  »t  encoun¬ 
ters  the  dense,  low-altitude  atmosphere.  The  third  example 
illustrates  .1  release  from  the  Space  Shuttle  orbiting  at  400-km 
altitude  with  an  eastward  velocity  of  7.7  km/s.  The  vapor  is 
cicvteJ  from  the  Shuttle  with  a  downward  vclocilv  of  .VI  km's 
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f'tg.  I I  I  rave!  line  .trnl  ih.tancc  for  refea.es  of  C  O.,  at  -toO-Mri  altitude  daring  solar  minimum.  Azimuths  of  0‘‘  and 
|M>'  ate  .tieii.i  and  against  the  mbit  velocity.  re.pecliv  eiy  Release.  »:ih  a  positive  elevation  angle  are  not  colli,  or.allv 
stopped  in  less  than  a  time  of  <»)  s  or  a  distance  of  t-  ■  P  fan. 


T  he  free-streaming  vapor  is  moving  at  it  high  velocity  (figure 
9c)  and  will  travel  Kir  many  tens  of  kilometers  before  being 
eolhsionally  stopped.  When  the  vapor  enters  a  diffusive  state, 
its  length  will  be  double  that  shown  in  the  figure. 

If  the  vapors  tire  released  from  rocket  engines,  it  is  unlikely 
that  only  one  species  will  he  released.  The  Space  Shuttle  en¬ 
gines,  for  instance,  release  substantial  quantities  of  CO,,.  N,. 
HjO.  and  If,.  The  fifth  example  deals  with  II,  vapor  expan¬ 
sion.  This  simulation  is  identical  with  the  fourth  simulation 
except  that  the  values  of  molecular  mass  and  collisional  cross 
section  for  H,  instead  of  CO,  are  used.  I  he  concentration  and 
velocity  for  this  release  are  illustrated  in  f  igure  10.  Note  that 
the  scale  in  figure  10  is  1.5  times  the  scale  in  figure  9. 

All  of  the  examples  dealt  with  the  release  of  lO-'*  molecules. 
1  o  find  the  results  for  the  release  of  different  quantities  of 
vapor,  scale  the  concentrations  according  to  the  amount  re¬ 
leased.  I  his  procedure  is  valid  as  long  as  the  release  quantity 
does  not  violate  the  criteria  established  by  f  igure  .1  in  section 


8.  Tkxvh  Dim  wc  i  ton  Scri  ksonk  Ru  t-  xsi  s 

As  demonstrated  in  the  previous  section,  supersonic  vapors 
may  come  to  rest  a  considerable  distance  Ironi  their  point  of 
relea.e  In  this  section,  the  travel  distance  for  the  bulk  of  the 
released  vapor  is  calculated  for  a  variety  of  injection  condi¬ 
tions. 

A  molecule  is  said  to  be  collisionally  halted  when  it  loses  a 
large  fraction  (say.  90'})  of  its  original  momentum  due  to 
molecular  collisions,  from  section  2.  the  mean  free  path  for 
supersonic  molecules  is  independent  of  velocity  hut  does  de¬ 
pend  on  the  concentration  of  the  background  atmosphere. 
Thus  the  condition  for  the  loss  of  90 of  the  momentum  of  a 
supersonic  molecule  is  approximated  by 


where  X  is  the  spatially  dependent  collision  frequency  and  yds 
is  the  integral  along  a  molecular  path.  In  a  uniform  atmo- 
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f  if)  13.  H,  travel  lime  and  distance  for  releases  under  the  same  conditions  as  I  igure  1 2  The  H,  molecules  are  coltish '-'alb 

halted  mote  rapidly  than  the  CO,  molecules. 


sphere,  a  typical  molecule  would  move  \  In  (10)  (2.3  mean  free 
path,). 

An  alternate  way  to  calculate  the  vapor  travel  distance  is  to 
consider  the  partitioning  of  the  molecules  between  the  natural, 
collisionless  response  [n^x,  /)].  and  the  forced,  collision-domi¬ 
nated  response  [rr.(x.  ()].  When  the  number  of  molecules  repre¬ 
sented  by  /7.(.v,  /)  drops  below  I0rc,  of  the  amount  released. 
9t.‘  ~t  of  the  pas  will  he  in  a  collision-dominated  state.  The  time, 
i;.  when  this  occurs  is  found  from 

J  n,(x.  r,)i/x3  =  QJ 10 

The  travel  distance  to  the  point  where  this  occurs  is  the  initial 
source  velocity  multiplied  by  neglecting  the  effects  of  grav¬ 
ity 

Doth  methods  yield  similar  results  for  the  vapor  travel  dis¬ 
tance  The  later  formulation  involving  the  spatial  integration 
o!  the  natural  response  is  used  here. 

The  following  calculations  consider  releases  of  1  x  1  (V* 
molecules  being  ejected  for  the  Space  Shuttle  moving  with  a 
horizontal  velocity  of  7.7  km/'s  The  injection  velocity  relative 
to  the  Shuttle  iv  .3.1  km  s.  The  travel  distance  is  shown  to  be  a 
I'.maion  of  the  background  atmospheric  concentration,  the 
injection  azimuth  and  elevation,  and  the  type  of  species  re¬ 
leased  In  the  release  coordinate  system,  the  elevation  is  mea¬ 
sured  from  the  horizontal  plane,  and  the  azimuth  is  the  angle 
between  the  Shuttle  velocity  vector  and  the  projection  in  the 
horizontal  plane  of  the  exhaust  injection  sector 

During  solar  minimum,  the  concentration  of  the  nevitr.it 
atmosphere  it  TOO  km  altitude  is  only  3.2-1  v  1  *  in  1  (see 

I  able  2).  Consequently,  high-velocity  sapor  releases  require 
the  maximum  time  and  travel  distance  to  he  colhsionally 
stopped  during  this  period  figure  I  I  illustrates  the  travel  lime 
.n. d  distance  for  the  sample  point  release.  Releases  with  a 
d  '  ■'  nw .ird  v  clouts  component  require  a  minimum  of  20  s  and 
2r«)  km  to  stop  This  corresponds  to  74  7  km  decrease  in 
altitude  and  1X5  A. km  movement  horizontally  tor  a  release  at 

Vi  elevation  |  he  travel  distance  for  releases  at  zero  degree 
elevation  is  close  to  (•< >0  km.  nearly  independent  of  azimuth 
I*  ; , .  *  rl  a.t'o;:  angle- v  ■■  a  I  ;  r  i .  a!  Pi  t  an.  os  er'*  'ter  t!  in  o'  ■ 


km  because  the  vapor  moves  into  more  tenuous  regions  of  the 
atmosphere  at  higher  altitudes. 

for  the  maximum  solar  cycle  conditions  given  in  Table  2. 
the  CO,  vapor  stop  times  and  distances  are  greatly  reduced 
(figure  12).  The  travel  lime  has  a  minimum  of  4.9  s  for 
releases  at  -30°  elevation  and  05  azimuth  The  minimum 
travel  distance  (47  km)  occurs  for  -50°  elevation  and  IkO* 
azimuth.  Supersonic  releases  rn  the  horizontal  plane  travel  a 
fixed  distance  (2.3  mean  free  paths)  before  they  are  colli- 
sionally  stopped.  Consequently,  at  zero  degree  elevation,  the 
trasel  time  lias  u  strong  azimuthal  component,  but  the  travel 
distance  is  nearly  constant. 

The  calculations  for  solar  maximum  atmosphere  are  re¬ 
peated  with  II,  instead  of  CO,  (figure  13).  The  light  H- 
molecule  rapidly  loses  momentum  to  the  heavier  background 
gas.  The  1 1,  sapor  is  stopped  w  iihin  25  km  ofihe  release  point, 
a  factor  of  2  less  than  r he  CO,  travel  distance. 

9.  StMMART  AND  Cost  I  l  MUSS 

The  theory  outlined  above  is  applicable  to  a  wide  range  of 
gas  release  experiments.  The  analysis  is  specifically  designed  to 
model  high-velocity  (i.e.,  Space  Shuttle)  releases  used  for  iono¬ 
spheric  depiction  experiments.  The  theory,  however,  tan  he 
applied  to  experiments  involving  the  injection  of  neutral  va¬ 
pors  for  use  as  tracers,  instability  triggers,  and  gravity  wave 
generators. 

The  theory  may  be  expanded  in  several  areas.  The  effects  of 
the  high-velocity  releases  on  the  background  atmosphere 
could  he  calculated  simultaneously  solving  the  coupled  Boltz¬ 
mann's  equations  (4)  and  (5)  Colli'ional  ionization  and  dis- 
assoeiation  may  be  included  in  the  modeling  by  adding  a 
velocity-dependent  loss  term  to  (6).  Chemical  losses,  due  to 
reactions  between  the  injected  gas  and  the  atmosphere,  may  be 
similarly  handled. 

I  he  present  analysis  has  predicted  several  fundamental  ef¬ 
fects  on  the  sapor  as  it  moves  supersonically  through  a 
tenuous  neutral  atmosphere  The  release  will  become  elon¬ 
gated  due  to  molecular  scattering.  During  the  collision;!!  proc¬ 
ess  heating  will  occur  The  m  .latitude  of  the  heating  depends 
mi  the  release  velocity  The  v  -par  travel  distance  is  a  i >« - ■ 


4  '^4 


Hi  KMUk:»i  Mk.ii  A i  mi  t»t  G  vs  K  1 1 1  xs i  s 


of  the  background  atmospheric  concentration  anti  of  the  direc¬ 
tion  id  (lie  vapui  injection. 
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D.  RADIOWAVE  PROPAGATION 


Propagation  of  radiowaves  in  the  ionosphere  can  be  substantially 
changed  due  to  the  injection  of  rocket  exhaust  vapors.  Artificially 
produced  plasma  holes  have  strong  plasma  gradients  which  affect  radio 
ray  trajectories.  We  have  examined  propagation  in  the  modified  iono¬ 
sphere  for  frequency  ranges  including  VLF  (3-30  kHz),  HF  (3-30  MHz), 
and  VHF  (30-300  MHz). 

/-t  VLF  waves  propagate  through  the  ionosphere  in  whistler  mode. 

Plasma  holes  will  defocus  these  rays  causing  a  reduction  in  signal 
strength.  As  plasma  flows  from  above  into  ionospheric  holes,  field 
aligned  reductions  in  plasma  concentration  are  formed.  These  depletion 
tubes  may  act  as  guides  for  VLF  waves  in  the  magnetosphere. 

At  HF  the  plasma  hole  focuses  transionospheric  wave  energy.  The 
magnitude  of  the  focusing  can  be  30  dB  or  more.  Ray  refraction  in  the 
bottom-side  ionosphere  is  changed  in  the  vicinity  of  a  plasma  hole. 
Depending  on  the  relative  locations  of  transmitter,  hole  and  receiver, 
the  changes  in  subionospheric  propagation  can  include  focusing,  defocusing, 
multipath  fading,  or  leakage  of  the  signal  out  of  the  hole. 

At  VHF  and  above  the  ionospheric  hole  produces  slight  focusing  of 
signals  propagating  from  satellites  to  the  ground.  Changes  in  phase 
and  polarization  of  VHF  waves  penetrating  the  hole  provide  experimental 
diagnostics . 

Little  is  known  about  scintillations  which  could  be  induced  by 
artificially  produced  holes  in  the  equatorial  ionosphere.  If,  as 
suggested  in  the  previous  section,  irregulari ties  form,  they  could 
introduce  random  phase  and  amplitude  changes  in  VHF  and  UHF  signals 
propagating  through  the  disturbed  region.  This  could  have  major  impact 
on  satellite  communications. 

The  changes  in  radiowave  propagation  are  described  in  greater  detail 
in  the  following  article  which  was  presented  at  the  AIAA  16th  Annual 
Meeting,  May  5-11,  1980,  in  Baltimore,  Maryland.  This  paper  has  been 
submitted  for  publication  to  the  Journal  of  Spacecraft  and  Rockets. 


CHANGE S  IN  THE  ELECTROMAGNETIC  PROf’E RT  I f  S  Or  Till: 
lll’I’LR  ATMOSPHERE  DUE  EO  ROCKET  EITLUKNTS 


Raul  A.  tier nhardt * 
Stanford  Electronics  Laboratories 
Stanford  University 
Stanford,  California 


Abstract 

Artificial  reduction  of  the  plasma  concentra¬ 
tions  in  the  ionosphere  and  magnetosphere  can 
affect  E-M  wave  propagation.  T fie  molecular 
species  found  in  rocket  exhaust  vapors  promote 
recombination  of  electrons  and  ions  (i.e.,  0*,  H*. 
and  He*)  in  the  upper  atmosphere.  This  recombina¬ 
tion  produces  localised  regions  with  strong 
variations  in  radio  refractive  index.  Regions  of 
ionospheric  plasma  reduction  produce  focusing  at 
HE  and  VIIE  frequencies  and  defocusing  at  VEf 
frequencies.  Subionospheric  propagation  in  the 
vicinity  of  an  artificial  plasma  hole  can  suffer 
the  effect'  of  multi -path  fading  and  defocusing  by 
reflection  at  the  sides  of  the  hole.  Should 
irregularities  form  because  of  the  plasma  gradients, 
amplitude  scintillations  will  be  introduced  in 
waves  passing  through  the  regie  it. 


II.  Molecular  Species  Present  in  Rocket  E xhajj s_t 

Exhaust  vapors  from  chemically  powered  rockets 
may  react  with  the  plasma  in  the  upper  atmosphere . 
This  factor  is  not  normally  considered  by  the 
designer  of  rocket  propulsion  systems.  The 
designer  considers  a  rocket  fuel  in  terms  of  per¬ 
formance  rather  than  in  terms  of  effect  on  the 
environment . 

figure  1  illustrates  a  comparison  of  pro¬ 
pellants  by  specific  impulse  ami  bulk  density. 

All  of  the  propellants  in  tins  figure  produce 
exhaust  species  which  reach  chemically  in  the 
ionosphere,  causing  a  reduction  in  plasma  concen¬ 
tration.  Designing  a  useful  chemical  propulsion 
system  which  does  not  eject  reactive  exhaust 
vapors  may  not  be  feasible. 


I  .  I  n  t  rodux  1 1  or 

The  earth's  plasma  atmosphere  can  tie  divided 
into  four  regions  with  approximate  Uiundame 
given  as  follows.  The  [)-region  is  at  altitudes 
between  70  and  40  km.  The  L -region  extends  from 
40  to  100  km  altitude.  The  F -region,  lying 
between  ldO  and  ROD  km,  contains  the  bulk  of  the 
ionospheric  plasma.  Above  the  iono'ghere  is  the 
magnetosphere  wh it h  extends  from  HUD  km  altitude 
to  over  10  earth  radii. 

The  firing  ot  rocket  engines  in  the  uppe> 
atmosphere  releases  vapors  which  can  temporarily 
neutralize  some  of  t  fie  plasma  in  the  F -region  and 
to  a  lesser  degree,  in  the  nagnetosph-re .  This 
prm  ess  produces  regions  of  localized  plasma 
reduction  called  pla'ma  "holes".  Kadi  i  waves 
propagat  mg  through  the  disturbed  regions  can  tie 
strongly  refracted  and  Mattered. 

The  obiect  ot  this  paper  is  to  survey  the 
variety  of  effects  that  artifii  l,il  1  y -produi  ed 
[i  I  a  sum  depletions  can  have  on  elei  tromagnei  it 
waves,  for  this  study,  numerital  models  ot  the 
upper  atmosphere  are  used  to  preditt  tfie  size  and 
durations  ot  the  depletons.  Iwo  and  three- 
dimensional  ravtraiing  programs  are  employed  to 
investigate  the  changes  in  the  radiowave  ray 
trajectories.  This  study  is  motivated  by  the  need 
to  evaluate  tFie  consegueru  es  of  temporary  alter, i- 
tiori  of  tfie  ionosphere .  Such  alteration  might  tie 
produced  by  the  launches  required  to  construct 
large,  orbiting  spare  Systems  such  as  the  Solar 
Rower  Satellite  (SRS). 

*  Research  Associate,  Kadiosc  lent  e  Laboratory 
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Tahir  I  lists  t ho  exhaust  s pec  ies  tor  two 
sol  id  rocket  propel  lant* .  Hu*  species  designated 
|»y  *  an'  known  to  produce  nod  1 1  i <  a t  1  on  ot  the 
urn  os  phere  .  T  hi  *  exhaust  t  run  solid  rocket  motors 
may  also  contain  chlorii.e  com;**  Hinds  whkli  tan 
effect  the  stratosphere  o/one  concentration. 

Most  solid  no* kets  burn  >nlv  in  the  lower  atmo¬ 
sphere  and  do  not  afteit  t  tie  mno sptiere . 


[  xhaust  vapors  from  liquid  rocket,  empties 
usual  I,  contain  tower  spot  ies  than  are  produced  by 
solid  rockets  (Tahh*  II).  This  exhaust  is  mainly 
composed  of  oxides  and  hydrogen  and  carbon  as  well 
as  unhurnrd  molecular  hydrogen.  Most  of  these 
gases  will  produce  plasma  depletions  if  injected 
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Table  I. 


Typical  Solid  Rocket  Exhaust 
Species  and  Concentrations 


Oxidizer 

hh4ciu4- 

•  -66 

Double-base 

Binder 

ich2: 

1  . 

-- 

Fuel 

ai 

Al 

-22. 

Combustion 

temperature  ( )  .>644 

4017 

Exhaust 

fi 

0.0233 

H 

0.1036 

products , 

Cl 

0 . 0084 

0 

0.0004 

fraction 

*011 

0 . 0006 

*011 

0.0046 

HC1 

0 . 54o? 

*h2 

0.7852 

*H.; 

1 .3183 

*ll?0 

0.1404 

*Ho0 

0.1730 

CO 

1.0910 

CO 

0.8339 

*C02 

0.0295 

*C0., 

0.0192 

*N0 

0.0004 

*N, 

0.2869 

*N2 

0.7951 

A1C1 

0.0011 

Al 

0.0002 

A1C1, 

0.0028 

A10 

0.0001 

AIC13 

0.0001 

ai2o 

0.0004 

A12°3 

0.4057 

A12°3 

0.4071 

*  Promotes 

ionospheric  plasma 

reduction . 

Table  II.  Typical  Liquid  Rocket  Exhaust 
Species  and  Concentrations 


Oxidi zer 

N2°4 

°2 

Fuel 

MMH 

h2 

Mixture  ratio  (0/F) 

1 .65 

6.0 

Theoretical  Is,  sec 

313 

388 

Combustion 

temperature  ("K) 

3250 

2869 

Exhaust  products,  *(!., 

0.2411 

0.2500 

mole  j,,.4' 

0  0.2738 

0.7500 

fraction  2 

CO 

0.0501 

0 

*ro 

2  0.1216 

0 

*N, 

0.3134 

0 

Promotes  ionospheric  plasma  reduction. 


into  the  f -region  ionosphere. 


The  spatial  extent  of  the  ionospheric  modifi¬ 
cation  is  affected  by  1)  the  total  amount  of 
exhaust  released  and  .')  the  manner  the  exhaust 
vapors  are  layed  downl.2.  The  amount  released 
depends  on  t  tie  eriqine  flow  rates  and  the  duration 


ot  the  enijine  burn.  The  flow  rate  of  the  Space 


-1 


Shuttle. Orbital  Maneuvering  Subsystem  is  8.71  kq S 
engine’  and  of  t.he  Space  Shuttle  Main  Enqine  (SSMl  ) 
T6H.4  kg  S’  enqine"  .  The  proposed  heavy  lift 


launch  vehicle  (IILLV)  for  the  SI’S  uses  14  SSME's, 
producin')  a  flow  of  65S7.6  kcj/sec. 


III.  Chemical^  Mei  han  mans  for 

’tonosjJherjc  Hole  formation 

A  complete  description  of  ionospheric  hole- 
formation  involves  chemistry,  gas  dynamics,  diffu¬ 
sion  and  thermal  processes.  The  details  of  these 
processes  can  be  found  in  several  reports-**"*. 

Briefly,  plasma  holes  are  formed  by  chemical 
reaction  of  vapors  in  rocket  plumes  with  ions  in 
the  F-layer.  Normally,  the  0*  j0n,  which  makes  up 
the  bulk  of  the  ionosphere,  is  removed  by  the 
following  reaction. 

O4  *  N2  -  N0+  +  N  .  rate:  1.3  x  10"13  cm3/ sec 

The  rate  of  this  reaction  is  at  least  1000  times 
lower  than  the  rates  for  the  reactions  between  0 
and  most  rocket  exhaust  species. 

Consider  competing  reactions  with  rocket 
exhaust 

0+  +  -*■  0H+  +  H  rate:  2  x  10 cm3/sec 

0+  +  H^O  H^0+  +  H  rate:  2.3  x  10  ^  cm3/sec 

0  +  C02  +  02  +  CO  rate:  1.2  x  10  ^  cm3/ sec 

The  relatively  large  rates  for  these  reactions 
cause  the  rapid  removal ^of  the  0+  ion.  The 
molecular  ions  (N0+,  OH  ,  H-0.  and  02+ )  quickly  - 
recombine  with  electrons  at  rates  around  2  x  10"/ 
cimVsec.  Similar  reactions  can  cause  reduction  of 
the  H+  and  He+  ions  in  the  magnetosphere. 

After  the  exhaust  vapors  are  released,  they, 
will  disperse  mainly  through  diffusive  expansion  . 
The  lifetime  of  the  ionospheric  hole  is  influenced 
by  transport  of  ions  into  the  affected  region  as 
well  as  by  re-ionization  of  the  neutrals  by  the 
sun’s  extreme  ultraviolet  radiation.  The  size  of 
the  hole  is  controlled  by  the  expansion  and  trans¬ 
port  of  the  injected  exhaust  vapors.  All  of  the 
processes  are  contained  in  the  numerical  models  we 
have  used  in  our  studies3'1*. 

The  models  are  used  to  generate  estimates  of 
disturbances  produced  by  the  release  of  vapors  in 
the  ionosphere.  The  changes  in  radio  wave  propaga¬ 
tion  are  investigated  by  tracing  rays  through  the 
perturbed  ionospheric  regions.  The  rest  of  this 
paper  describes  typical  effects  of  artificial 
ionospheric  holes  on  radio  waves  propagating  at  VLT 
(3-30  kHz),  at  HF ( 3-30  MHz,  at  VIIF(30-300  MHz)  and 
above . 

]_V_.  Modification  of  Propagation.  Vlf 
Waves  in  tin’  Ionosphere 

Long  distance  propagation  of  very  low  fre¬ 
quency  (VLT)  waves  may  be  either  I)  via  the  earth- 
ionosphere  waveguide  or  2)  via  field  aligned  ducts 
in  the  earth's  magnetosi  here .  The  first  type  of 
propagation  employs  trie  earth’s  surface  and  the 
lower  edge  of  the  ionosphere  for  guiding  the  waves. 
Anv  modification  of  the  ionosphere's  lower  boundary 
will  change  the  characteristics  of  the  earth-iono¬ 
sphere  waveguide. 

The  lower  ionosphere  is  composed  of  positive 
polyatomic  ions  such  as  NO4,  ().,*  (in  the  D  and  E 


region),  negative  polyatomic  ions  such  as  '+(D 
region)  and  hydrated  positive  ions  such  as  H-H,0 
{0  region).  These  ions  are  not  nearly  as  sensitive 
to  the  chemical  changes  as  are  the  0+  ions  of  the 
F-layer,  Consequently,  the  effect  of  rocket 
exhaust  on  propagation  via  the  earth-ionosphere 
waveguide  is  expected  to  be  small. 

The  second  type  of  VLF  propagation  (called 
whistler  mode)  is  illustrated  in  Figure  2.  A  VLF 
electromagnetic  wave  is  transmitted  from  a  source 
into  the  earth-ionosphere  waveguide.  Some  of  the 
wave  energy  leaks  upward  through  the  bottomside 
ionosphere,  traveling  in  an  unguided  mode.  If 
this  energy  is  captured  by  an  enhanced  column  of 
plasma  (called  an  enhancemend  duct),  the  wave  can 
be  guided  along  the  earth's  magnetic  field  lines 
to  the  conjugate  hemisphere.  Amplification  may 
occur  as  the  wave  interacts  with  the  radiation 
belt  particles  in  the  magnetosphere^ _  Radio  rays 
leaving  the  duct  will  travel  in  an  unducted  mode 
to  the  bottomside  ionosphere.  If  the  ray  direction 
is  nearly  perpendicular  to  the  ionosphere,  signals 
will  couple  into  the  earth-ionosphere  waveguide 
and  may  be  received  on  the  ground.  Further  details 
on  whistler  mode  propagation  can  be  found  in  the 
text  by  Hell  iwel  P. 
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Fig.  2  Schematic  of  the  magnetospheric 

propagation  circuit  for  VLF  waves. 

Modification  of  the  F-region  ionosphere  will 
change  the  coupling  between  the  earth-ionosphere 
waveguide  and  the  ducts  in  the  magnetosphere.  This 
is  demonstrated  by  tracing  5  kHz  rays  through  an 
ionosphere  that  was  modified  by  the  release  of  100 
kg  of  H,  at  300  km  altitude.  The  electron  density 
contours  one  minute  after  the  H-  release,  are 
illustrated  in  Figure  3.  i 

For  rays  to  be  captured  by  a  whistler  duct, 
their  wave  normal  components  must  be  nearly 
aligned  with  the  earth's  magnetic  field.  At  high 
geomagnetic  latitudes,  the  magnetic  field  lines 
are  nearly  vertical.  Upward  propagating  VLF  waves 
have  their  wavenormals  refracted  vertically  as  they 
enter  the  ionosphere  at  the  lower  boundary.  As 
these  waves  propagate  in  a  horizontally  stratified 
ionosphere,  the  wavenormals  tend  to  remain  vertical 
and  the  waves  tend  to  couple  into  any  enhancement 
ducts  they  encounter. 


NIGHT-TIME  ELECTRON  CONCENTRATION 
LOG.0(Number  Density, cm  ') 
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Fig.  3  Ionospheric  hole  produced  by  the  point 
release  of  100  kg  of  at  300  km 
altitude.  One  minute  after  the  release, 
the  central  plasma  concentration  is 
reduced  by  a  factor  of  10. 

The  wavenormals  of  the  rays  are  sharply  bent 
by  the  plasma  gradients  of  an  ionospheric  hole. 
Figure  4  illustrates  the  bending  of  b  kHz  rays  by 
the  ionospheric  hole  of  Figure  3.  The  rays  are 
emitted  from  a  ground  station  located  at  5b°  geo¬ 
magnetic  latitude.  The  three  rays  to  the  left  of 
the  hole  are  unaffected.  In  the  vicinity  of  the 
hole,  the  rays  are  defocused  and  their  wavenormals 
are  bent.  (The  wavenormal  directions  are  indicated 
by  the  arrows  at  the  end  of  the  ray  trajectories). 
These  rays  cannot  couple  into  ducts  in  the  magneto¬ 
sphere. 
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Fig.  4  Defocusing  of  5  kHz  rays  by  an  iono¬ 
spheric  hole  centered  on  the  small 
c i rc  1  e . 
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A  reciprocal  effect  occurs  for  signals  exit¬ 
ing  a  duct  over  an  ionospheric  hole.  Figure  5 
illustrates  the  defocusing  and  wavenormal  bending 
for  5  kHz  rays,  propagating  downward.  Without  the 
hole,  all  rays  would  cross  the  lower  ionosphere 
and  propagate  to  the  earth's  surface.  Between 
53. 6°  and  Bo5  geomagnetic  latitude,  the  iono¬ 
spheric  hole  prohibits  transmission  to  the  ground 
of  the  5  kHz  wave. 


Fig.  5  Refraction  of  downward  proiagating 
rays  by  an  F-region  depletion.  A 
transmission  coefficient  of  zero 
indicates  rays  which  do  not  penetrate 
the  bottoms ide  ionosphere. 

Figure  6  illustrates  the  required  wavenormal 
angle(s)at  500  km  altitude  for  propagation  through 
the  ionospheric  hole  to  the  earth's  surface.  The 
wavenormal  direction  is  defined  with  respect  to 
the  vertical.  The  dashed  portion  of  the  curve 
shows  where  strong  defocusing  occurs.  Rays  leav¬ 
ing  enhancement  ducts  have  wavenormal  angles  usu¬ 
ally  less  than  15°.  In  regions  where  the  pene¬ 
tration  wavenormal  angle  is  greater  than  this 
value,  coupling  between  whistler  ducts  and  the 
earth-ionosphere  waveguide  does  not  occur. 

Fading  of  the  signal  amplitude  will  occur  if 
two  or  more  rays  reach  the  ground  receiver,  each 
ray  having  propagated  on  a  disjoint  path.  Between 
53. 6°  and  53.9°  latitude,  rays  with  three  differ¬ 
ent  wavenormal  angles  can  reach  the  earth's  sur¬ 
face.  Consider  a  VIF  satellite  transmitting  a 
5  kHz  omnidirectional  signal  from  a  500  km  orbit. 
When  this  satellite  enters  the  fading  region 
(Figure  6),  amplitude  fluctuations  will  be  per¬ 
ceived  on  the  ground. 


FADING 


Fig.  6  Required  wavenormal  angle  at  500  km 
altitude  for  ground  reception  of  a 
5  kHz  wave.  The  angles  are  defined 
with  respect  to  the  vertical. 

One  further  way  that  VLF  propagation  in  the 
magnetosphere  may  be  modified  by  exhaust  vapors  is 
with  the  creation  of  artificial  ducts.  When  the  F- 
layer  is  locally  depleted,  plasma  flows  along  mag¬ 
netic  field  lines  from  the  magnetosphere  into  the 
depleted  region.  This  flow  will  produce  a  tube  of 
reduced  plasma  concentration  extending  up  into  the 
magnetosphere.  This  "depletion  duct"  takes  three 
or  more  hours  to  form. 

Bernhardt3  has  made  estimates  of  the  amount  of 
reduction  in  the  content  of  a  magnetic  flux  one 
wouid  expect  from  vapor  releases.  For  instance, 
1000  kg  of  H_  released  at  300  km  altitude  into  the 
daytime  ionosphere  is  calculated  to  produce  10o 
reduction  in  the  electron  content  of  a  tube  cen¬ 
tered  along  a  L=4  (earth  radii)  magnetic  field  line. 
The  radius  of  this  tube  will  be  650  km  at  the 
equator. 

This  type  of  duct  may  be  able  to  guide  VLF 
waves  in  the  magnetosphere.  As  mentioned  pre¬ 
viously,  enhancement  ducts  are  efficient  guides  of 
VLF  waves.  Depletion  ducts  can  also  guide  VLF, 
using  one  of  two  modes.  The  first  mode  involves 
large  wavenormal  angles  for  rays  which  do  not 
penetrate  the  lower  ionosphere  to  the  ground.  The 
second  mode,  invol ves, trapping  at  the  lower  edge 
of  the  depletion  duct  .  This  "one  sided"  mode  can 
be  excited  from  the  ground  for  guiding  a  wave  to 
the  conjugate  hemisphere.  A  depletion  duct  of  10i- 
should  be  sufficient  for  this  type  of  guiding. 

V.  Artificial  Changes  in  High 
Frequency  tT-DO  MHz')'  Waves. 

The  peak  concentration  of  the  F-layer  defines 
a  lower  limit  for  the  frequency  of  penetration  of 
HF  waves.  This  frequency  is  given  by  f  F?  = 


\/'> 

[80.6  N  1  Hz  where  N  , 

L  max  -  max 

tration  in  electrons  m  .  An  ionospheric  hole. 


is  the  peak  concen- 


located  at  the  F-layer  peak,  will  decrease  f  r,’, 
allowing  lower  frequency  (HF)  waves  to  penetrate. 


This  effect  was  proposed  by  Papagiannis  and 
Mendillcr  for  use  in  low  frequency  radio  astronomy. 


Refraction  of  an  HF  wave  by  the  plasma 
gradients  of  the  ionospheric  hole  will  greatly 
distort  the  propagating  phase  front.  The  magni¬ 
tude  of  this  distortion  is  a  function  of  the  wave 
frequency,  the  plasma  gradients  in  the  hole  and 
the  direction  of  wave  propagation  through  the  hole. 


Figure  7  is  an  example  of  downward  propaga¬ 
ting  HF  rays  throuqh  the  ionospheric  hole  shown 
in  Figure  2.  The  rays  are  superimposed  on  con¬ 
tours  of  constant  electron  concentration.  The 
ordinary  left-hand  circularly  polarized)  mode  is 
designated  0-mode.  The  effect  of  the  reduction  in 
plasma  concentration  is  to  focus  the  downgoing 
waves.  This  is  in  contrast  with  the  defocusing  of 
VLF  waves  discussed  in  the  previous  section.  In 
the  example  shown  in  Figure  7,  the  focal  point  for 
the  7  MHz  0-mode  lies  on  the  ground.  The  iono¬ 
spheric  "lens"  produced  by  a  100  kg  release  of 
hydrogen  at  300  km  altitude  is  roughly  equivalent 
to  a  perfect  lens  with  a  diameter  of  36  km.  The 
use  of  such  a  lens  for  radio  astronomical  observa¬ 
tions  has  been  proposed  by  Bernhardt  and  da  Rosa°, 
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Fig.  7 


Propagation  of  a  7  MHz  ordinary  wave, 
throuqh  the  ionospheric  hole  illus¬ 
trated  in  Figure  3.  The  downward 
propagating  rays  are  focused  at  the 
enc i rc led  point. 


Ionospheric  plasma  holes  can  also  affect  the 
propagation  of  high  frequency  waves  between  two 
points  on  the  earth's  surface  (called  subionn- 
spheric  propagation) .  Depending  on  the  geometry 
of  the  transmitter,  receiver  and  ionospheric  hole, 
the  HF  signals  can  experience  focusing,  defocusing 
and/or  fading. 


Representative  effects  of  ionospheric  holes 
on  HF  propagation  are  taken  from  the  work  of  Helms 
and  Thompson^.  They  traced  rays  through  an  iono¬ 
spheric  perturbation  which  we  estimate  could  be 
produced  by  the  release  of  3000  kg  of  H,  into  the 
nighttime  F-layer.  In  their  work,  t tie  electron 
concentration  is  normalized  by  square,of  the 
frequency,  as  follows:  X  -  80.6  ’0^/ ic  where  f  is 
the  wave  frequency  Hz  and  Ng  is  the  electron 
concentration  in  i'i  .  Xmax  is  the  value-at  the 
F-layer  peak  given  by  Xl;;ax  =  80.6  Nmax/i  • 

The  following  is  a  sutnmary„of  the  discussion 
presented  by  Helms  and  Thompson..,  With  -  4 
(for  example  with  Nn,i1x  -  1  x  10  1  m'3  and  1  = 

4.5  MHz),  the  unperturoed  ionosphere  produces  the 
0-mode  trajectories  shown  in  Figure  8.  The  iono¬ 
spheric  holes  modify  the  trajectories  as  shown  in 
Figure  9.  The  center  of  the  hole  is  indicated  by 
a  dashed,  vertical  line. 
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Fig.  8  Ray  path  trajectories  in  an  undis¬ 
turbed  ionosphere  with  Xmax  =  4  (from 
Helms  and  Thompson 3). 

Without  the  hole,  the  point  of  maximum  ray 
altitude  tor  rays  with  20"  to  40"  elevation  angles 
is  230  km  from  tne  transmitter.  With  the  hole 
placed  440  km  away  from  the  transmitter  (Figure  9a) 
the  same  rays  propagate  greater  distances  reaching 
maximum  altitudes  as  much  as  600  km  away  from  the 
transmitter.  This  could  have  a  desirable  effect 
cn  long  distance  communications. 

Moving  the  hole  closer  to  the  transmitter 
affects  rays  with  higher  intial  elevation  angles. 
Tor  hole-transmitter  separations  less  than  333  km, 
some  of  the  rays  are  reflected  back  to  the  trans¬ 
mitter  (Figures  9b  and  c).  finally,  for  separa¬ 
tions  less  than  120  km,  some  rays  escape  through 
the  ionospheric  hole  (Figure  9c).  When  the  hole 
is  directly  over  the  transmitter  (Figure  9dl  the 
rays  are  reflected  once  or  twice  try  the  sides  of 
the  hole  before  reaching  the  ground. 

In  some  cases,  an  iono  oheric  hole  can  pro¬ 
duce  focusing  or  defocusing  of  HF  waves.  Reflec¬ 
tion  from  the  side  of  the  hole  tends  to  spread 
the  ray  trajectories  (Figure  10a).  By  lowering 
the  wave  frequency  from  4.5  MHz  (Xmax  ~  4)  to  2.8 
MGz  (*max  =  1°)  focusing  of  the 'reflected  rays 
can  be  produced  (Figure  10b). 
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Fig.  9  Ray  path  trajectories  as  a  function  of  location  for  a  hole  with  a  depression  ration  of 
0.1,  a  halfwidth  of  22 2  km,  and  Xinax  =  4  (From  Helms  and  Thompson9). 


(a)  Fiefocusing 


Focusing  (b) 


Fig.  10  Focusing  and  defocusinq  effects  due  to  propagation  through  an  ionospheric  hole 
(From Helms  and  Thompson9). 


An  attempt  was  made  to  observe  perturbations 
ot  HF  communications  by  rocket-induced  ionospheric 
holes.  On  Sept.  20,  1979,  a  large  hole  was  created 
by  the  Centaur  stage  of  a  rocket  placing  the  HEA0- 
C  satellite  into  orbit  (see  Mendillo'^).  Radio 
amateurs  were  solicited  to  participate  in  an  ex- 
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periment  to  observe  the  effects  of  the  hole  on 
communications  .  Beacon  transmissions  3.6,  7.1, 
14.1,  21.1,  2B.1  and  50.1  MHz  were  made  from 
Puerto  Rico.  Signal  strength  measurements  from 
over  150  amateurs  in  the  continental  United  States 
are  reported.  The  reports  indicated  that  some 


'j 


sum !  1  coi’iminica tii'ii  effects  may  have  been  produced. 
These  effects,  however,  were  of  the  some  magni  tude 
that  one  wojlJ  expert  with  the  ionosphere's  norma  I 
variahi  1  i  tv.  the  co-nun  i  cat  ion  effects  of  the 
Ht-V'-C  launch  are  still  under  study. 
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Fro-;  the  ‘heoreticjl  examples,  we  see  that 
spheric  holes  can  produce  a  variety  of  changes 
.,C'f  anospneric  Hi  communications.  Holes 
j r  1  >  aver  a  transmitter  will  defocus  the  wave 
with  large  elevation  ancles.  Under 
special  iced  conditions,  holes  several  nun- 
ore  ter  s  fro:::  a  transmi  tter  can  produce 
resultin'  in  an  enhancement  of  siunal 
e.  Multi -path  propagation  produced  by  the 
cause  fadinj. 
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Hi  the  wave  frequency  is  increased  above  JO 
‘‘Hz,  the  refraction  by  the  ionospheric  gradients 
become ?  less  pronounced.  Above  100  MHz,  the  ray 
bond  in:  is  so  weak  that  propagation  between  two 
points  on  the  ground  (via  ionospheric  reflection) 
becomes  unlikely.  Consequently,  the  only  effects 
of  interest  occur  tor  trvis-ionospheric  propaqation 
An  example  of  this  mi  :ht  be  si  ana Is  propagating 
tivi  a  satell ite  through  the  ionosphere  to  the 
'jr.vis.l. 


g  u 

Fiq.  11  Electron  .on tent  c nances  associated 
with  trie  launch  of  HE.AA-C.  The 
launch  toon  place  at  Deed  JT.  T no 
rocket  exhaust  first  crossed  the  .of 
propaqation  path  ten-ina t i r.q  on  tie 
Crooked  Stiver  station.  The  measure¬ 
ments  were  cade  us  in.'  a  VnF  signal 
transmitted  from  tne  STRIP  satellite 
located  in  .eosynerrmnunous  orbit  -it 
14.77  w  lc,ngi  tude,  7.J7  S  latitude. 


Focus  in,;  will  occur  for  VHF  (30  -  300  MHz) 
si  fils  travelling  vertically  through  ar,  ionospheric 
hole.  The  phenomenon  is  similar  to  that  shown  in 
Figure  7  but  to  a  much  lesser  degree.  For  fre¬ 
quencies  above  30  MHz,  the  amplitude  increase  on 
the  jround  is  approximately  given  oy 

fi  1  +  J~z 
l  -  a 

where  -  is  the  wavelength,  z  is  the  altitude  of  the 
hole,  a  is  tne  halfwidth  of  the  hole,  ,  is  the 
phase  charue  through  the  center  of  the  hole  and  (! 
is  the  aripl  i  tude, gain  resulting  from  propaqation 
through  the  hole  ’.  Ttie  phase  perturbation  can  be 
found  from  the  integrated  vertical  electron  content 
with  and  without  the  hole  , 

=  M.4..x  JP-.-JJ 

whore  71  is  the  reduction  in  electron  content  (in 
nr-)  through  tne  center  of  the  hole  and  f  is  the 
wave  frequency  (in  Hz). 

As  an  example,  consider  an  ionospheric  hole 
centered  at  an  altitude  of  300  km  with  a  content 
reduction  7,1  -  6  x  lo'b  in'"  and  a  size  a  -  20  km. 
For  a  wave  frequency  f  =  140  MHz,  one  obtains  an 
amplitude  gain  r,  -■  1.1.  This  .8  dB  gain  is  small 
compared  with  naturally  occurring  variations  in 
signal  amplitude. 

The  reduction  in  electron  content  associated 
with  ionospheric  holes  produces  changes  in  the 
prrv, i'  and  polarization  of  VHF  waves.  Measurements 
o'"  these  changes  have  been  made  by  a  number  of 
worker^  to  dej.5et.jind  study  artificially  produced 
hole-,'’  •  1  ’  J.  Recent  measurements  made 

after  the  launch  of  the  IILA0-C  satellite  are  shown 
in  Figure  II.  IHo  data  for  this  figure  was  pro¬ 
vided  by  Stanford  University,  Boston  University, 

Air  Form  Geophysical  Laboratory,  and  the  Depart¬ 
ment  of  the  Army,  Fort  Monmouth.  More  details  can 
be  found  in  the  paper  in  this  session  by  Mendillo'P 


The  previous  discussion  con-.er, 'rated  on-' 
ations  in  radiowave  propagation  due  to  tne  t~. 
plasma  gradients  which  for"  tne  plasma  hole, 
scale  size  of  these  gradients  is  the  order  of 
of  kilometers  or  more.  One  should  consider, 
ever,  the  possibility  of  the  hole  breaking  j; 
irregylari t ies.  As  suggested  b,  Anderson  and 
hardt"  and  by  O.sakow,  et  al.'-,  holes  create 
the  Potto-  side  T-region  at  the  .magnetic  equal 
become  unstable.  This  phenomena  ocu-r  natur 
ygd  »:as  ij’c’cn  reported  by  '  rubber  of  re.,e-ar<jie 


O  t  T 
Tne 


L-ern- 
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The  instability  growth  rate  increases  as  the 
hot  toms ide  plavva  gradients  become  steeper.  This 
growth  rate,  however,  is  decreased  by  enhanced 
electron-ion  recombination  in  the  region.  In  the 
case  of  exhaust  rc-’eases,  the  mechanism  which 
steepens  the  bottomside  F-layer  also  increases  tne 
recombination  in  the  region.  Consequently,  the 
exhaust  vapors  may  damp-out  any  instabilities 
which  would  otherwise  form.  If  the  plasma  hole  and 
the  exhaust  vapors  become  separated,  the  instabil¬ 
ity  can  grow.  Electric  fields  and/or  neutral  winds 
may  be  adequate  for  transporting  the  hole  and  the 
exhaust  vapors  away  from  each  other. 

The  effects  of  the  irregularities  on  V HI  (and 
higher  frequency)  propagation  depends  on  the  size 
and  distribution  of  the  irregularities.  Naturally 
occurring  irregularities  have  been  known  to  pro¬ 
duce  severe  amplitude  and  phase  scintillations  in 
trans-ionospheric  oropagating  signals. 

Vi  I.  Concluding  Remarks 

Current  plans  for  the  Solar  Power  Satellite1'' 
call  for  the  multiple  launches  of  chemical  rockets 
for  transportation  of  cargo  and  crews.  The  Heavy 
Lift  Launch  Vehicle  ( HLLV ) ,  the  Personnel  Launch 
Vehicle  (I’LV),  and  the  Personnel  Orbit  Transfer 
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Venw'e  ,;vr,  '•  a.  .Vj'iiMt  a  .til’s  tan  ti  a  I  amount  of 
i-xnu.t  it,  ’is  ui  Mu  upper  atmosphere.  Themaqni- 
!...!<•  a'-!  t  's’  Jurat  i. 'ii  ot  thn  ionospheric  plasma 
raj...  ti  .r.'.'.iMsi  a.  the  urvjine  burns  depends  on 
M:e  ,  ,i r  r  i v  ii i  ; i  b..ru  prut  lies.  These  plasma  reduc- 
ti’iis  •  V  a'a’e  :  „ith  a  Hamper  ot  existin'] 

i  s’. \c !-  ‘  •  ••'■•-- I  .  The  re  true’  ion  effects  tv,  the 
t  1 : .  a  1  !  v  -e-s-Jui  ed  pi  J'ma-qradi ent s  can  be  reud- 
tl.  at.  ’  :  lis'-'d  ceupl  in-:  rav  trap  in.)  to  the 
ion."  p-erip  -ode!  . 

T'-e  -a;  '<•  a-  -'''taint,  in  tills  field  of  re- 
,ea-  is  *,M’-ber  or  not  exhaust  releases  mil 
tri-;  ;er  i  rre  ;jl  ar  1 1  ies  in  t'ie  equatorial  ionosphere. 
This  is  of  special  mter-.t  to  the  SI’S  program  be¬ 
cause  the  HLL %  are  expected  to  conduct  orbit 
s  i  iva  1  an  rat  ion  burns  at  the  equator.  The  amount 
of  material  released.  during  one  of  these  burns  may 
be  as  r.izn  as  9  x  1  O’’  kg  over  a  region  1 JOC  km 
luiid.  Should  tiie  HLIV  rake  an  expected  391  flights 
pc  year,  tiie  exhaust  vapors*-].'!.,'  continually  stimu¬ 
late  ionospheric  irregularities  which  could 
adverse!.,  if  feet  t  ra  ns  -  i  ono  spheric  propagation  in 
the  disturbed  region.  More  theoretical  and  experi¬ 
mental  work  needs  to  be  done  in  this  area. 

(Ac k r,p«l edne-ent '  This  work  was  supported  by  the 
Office  "ef  'Ja'vaT  Ses^ircii  Grant  NOOOl-T  -  77-C-llo5u 
and  Ar  uimie  National  laboratory  Contract  .31-109- 
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E.  ANALYSIS  OF  VLF  PHASE  AND  AMPLITUDE  DATA 


We  have  made  a  systematic  search  of  VLF  phase  and  amplitude  data  for 
rocket  launch  effects.  The  data  set  is  based  on  signals  received  at 
Byrd  Station,  Antarctica  after  having  propagated  through  the  earth- 
ionosphere  waveguide  from  one  of  several  Standard  Frequency  and  Time 
Stations  operating  at  VLF.  The  data,  taken  during  the  time  period  of 
September  1964  and  April  1970,  has  not  been  previously  analyzed.  Signals 
propagating  along  great-circle-paths  in  the  vicinity  of  Cape  Canaveral  , 
Florida,  were  examined  during  times  of  rocket  launches.  There  were  no 
systematic  fluctuations  in  the  phase  of  amplitude  that  could  be  associated 
with  rocket  launches.  From  this,  we  conclude  that  rocket  exhaust  releases 
in  the  lower  ionosphere  are  not  sufficient  to  produce  significant  fluctua¬ 
tions  in  VLF  signal  propagation  on  subionospheric  paths. 

We  did,  however,  isolate  periodic  fluctuations  in  the  data  that  were 
attributed  to  the  lunar  tide.  Also,  the  VLF  amplitude  data  was  shown  to 
contain  a  significant  increase  in  signal  strength.  The  effect  may  be  due 
to  power-line-induced  precipitation  of  energetic  electrons  which  ionize 
the  D-region  ionosphere.  The  details  of  this  analysis  are  given  in  the 
following  article  which  is  to  be  submitted  to  the  Journal  of  Geophysical 
Research . 


PERIODIC  FLUCTUATIONS  IN  THE  EARTH-IONOSPHERE  WAVEGUIDE 


Paul  A.  Bernhardt 
Kent  M.  Price 

Radioscience  Laboratory,  Stanford  Electronics  Laboratories, 

Stanford  University,  Stanford,  California  94305 

James  H.  Crary 

Jefferson  County  Data  Processing  Department 
Golden,  Colorado  90419 

ABSTRACT 

VLF  phase  and  amplitude  measurements  made  in  1967  at  Byrd  Station, 
Antarctica,  using  18.6  kHz  transmissions  from  NLK,  Washington,  have  been 
analyzed  for  lunar-induced  and  hebdomadal  variations.  The  VLF  propagation 
path  was  14,250  km  in  length,  and  nearly  constant  in  longitude.  The 
lunar-induced  variations  are  found  to  be  rlo:  in  phase  and  .-0.26  dB  in 
amplitude  and  are  explained  by  height  changes  in  the  D-region  of  the 
ionosphere.  The  hebdomadal  variation  in  VLF  amplitude  show-  a  0.46  dB 
increase  during  the  Sunday  daylight  hours,  relative  to  the  Monday-Saturday 
average,  and  is  ascribed  to  man-made  changes  in  the  D-region. 

INTRODUCTION 

The  phase  and  amplitude  of  VLF  signals  propagating  in  the  earth- 
ionosphere  waveguide  are  affected  by  fluctuations  in  the  D-region.  These 
fluctuations  can  be  either  irregular  transient  or  periodic.  Irregular  or 
transient  disturbances  have  been  associated  with  solar  flares,  magnetic 
activity,  meteor  showers,  whistler-induced  electron  precipitation  and 
geomagnetic  storms  (Chilton  et  al .  ,  1  964  ;  Potemra  and  Rosenberg,  1  973  ; 
Muraoka,  1979;  and  references  cited  in  these  works).  Sources  of  periodic 
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fluctuations  may  include  the  earth's  rotation  relative  to  the  sun,  the 
solar  27  day  rotation,  the  lunar  semi-diurnal  tide  and  the  hebdomadal 
variations  in  power  line  radiation.  In  this  paper,  only  the  periodic 
changes  in  VLF  amplitude  and  phase  are  considered. 

The  lunar-induced  semi-diurnal  variation  in  VLF  phase  has  been 
previously  investigated  by  Brady  and  Crombie  (1963).  They  found  that  the 
lunar  tide  produces  about  .*1°  of  variation  in  the  phase  measured  along 
propagation  paths  4000  to  5000  km  in  length.  This  phase  change  was  shown 
to  be  equivalent  to  a  *0.1  km  fluctuation  in  the  height  of  the  D-region. 
The  peak  height  of  the  D-region  was  found  to  occur  one  hour  before  upper 
and  lower  transit  of  the  moon.  It  had  been  previously  shown  by  Appleton 
and  Beynon  (1949)  and  by  Mi tra  (1955)  that  maximum  absorption  occurs 
approximately  one  or  two  hours  before  lunar  transit.  If  the  D-region 
produces  the  absorption,  it  seems  puzzling  that  the  peak  absorption  should 
occur  when  the  layer  is  elevated  to  its  highest  extent  by  the  lunar 
atmospheric  tide.  In  this  paper,  both  VLF  amplitude  and  phase  data  are 
analyzed  to  provide  further  information  about  the  lunar  tidal  effects 
in  the  D-region. 

Another,  more  recent  puzzle  concerns  the  effects  of  power  line 
radiation  (PLR)  on  the  energetic  particle  population  in  the  magneto¬ 
sphere.  A  number  of  papers  (Helliwell  et  al . ,  1  975;  Bullough  et  al  . , 

1976;  Park  1977;  Luette  et  al.,  1977;  Park  and  Helliwell,  1978;  Luette 
et  al . ,  1979;  Park  and  Miller,  1979;  and  Luette  et  al.,  1980)  have  sug¬ 
gested  that  PLR  may  1)  trigger  VLF  emissions  in  the  magnetosphere  and 
2)  change  the  trajectories  of  energetic  electrons  causing  particle 
precipitation  into  the  upper  atmosphere.  If  the  energetic  electrons 


2 


ionize  the  D-region,  they  might  produce  a  measurable  change  in  VLF  signals 
propagating  in  the  earth-ionosphere  waveguide.  Fraser-Smith  (1979)  and 
Park  and  Miller  (1979)  have  suggested  that  effects  due  to  PLR  may  be 
identified  by  looking  for  hebdomadal  periodicities  in  the  data.  The 
lower  usage  of  electric  power  on  weekends  (especially  on  Sunday)  could 
produce  variations  with  a  seven  day  period.  A  search  for  such  periods  is 
described  in  the  second  half  of  this  paper. 

VLF  AMPLITUDE  AND  PHASE  DATA 

The  paper  reports  on  the  analysis  of  VLF  amplitude  and  phase  measure¬ 
ments  made  in  1967  using  the  propagation  path  between  the  13.6  kHz  NLK 
transmitter  located  at  Jim  Creek,  Washington  (48.204°N,  -121.25°E),  and 
a  receiver  located  at  Byrd  Station,  Antarctica  (-79.97"N,  -120. 0°E).  The 
great  circle  propagation  path  is  shown  in  Figure  1.  The  length  of  the 
short  path  is  14,250  km,  and  the  length  of  the  long  path  (on  the  opposite 
side  of  the  Earth)  is  25,770  km.  The  short  path  has  a  nearly  constant 
longitude. 

Representative  plots  of  the  amplitude  data  are  illustrated  in  Figure 
2.  During  winter  (Figure  2a)  and  summer  (Figure  2c),  the  dawn-dusk 
terminator  crosses  the  VLF  propagation  path  at  a  nonzero  angle.  This 
produces  a  more  gradual  transition  between  day  and  night  values  than  is 
observed  during  the  periods  around  Equinox  (Figures  2b  and  2d). 

For  the  first  50  days  of  phase  measurements,  the  receiver  frequency 
was  6  x  10  ^  Hz  offset  from  the  NLK  transmitter  frequency.  This  produced 
a  gradual  shift  in  the  phase  measurements  (Figure  3).  After  the  receiver 
frequency  was  corrected,  the  phase  contained  relatively  slow  variations 
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on  which  the  daily  changes  are  superimposed.  Before  analysis  of  the  data, 
the  slow  variations  were  estimated  and  subtracted  from  the  phase  data. 

Like  the  amplitude  data,  the  winter  and  summer  phase  data  (Figures  4a  ana 
4c)  showed  more  gradual  transitions  between  day  and  night  tnan  aid  the 
measurements  taken  during  the  periods  around  Equinox  (figures  4b  and  4dJ. 

MODtL  OF  THE  EARTH- IUNUSPHLRL  ( L - 1 )  WAVEGUIDE 

The  E-I  waveguide  consists  of  a  spherical  lower  bounder,'  (i.e., 
the  earth's  surface)  and  an  upper  boundary  (i.e.,  the  D-region  at  the 
bottom  of  the  ionosphere).  Since  most  of  the  propagation  path  is  over 
the  ocean,  the  lower  surface  is  taken  to  have  infinite  conductivity.  The 
D-region  conductivity  parameter,  as  defined  by  Wait  and  Spies  (1%4;,  is 
given  by 

2 

ur  =  ;°  =  (2.5  ■  io5  s  ])  exp[f(z  -  h)]  ,  (1) 

where  u>  is  the  plasma  frequency,  v  is  the  effective  collision  frequency, 

3  is  the  bottom-side  normalized  conductivity  gradient,  h  is  the  reference 
height,  and  z  is  the  altitude.  Typical  coefficient  values  for  the  day¬ 
time  D-layer  are  B  =  0.3  km  ^  and  h  -  70  km.  The  conductivity,  .  ,  is 

■) 

directly  proportional  to  electron  concentration  because  of  the  / 
dependence  and  inversely  proportional  to  neutral  concentration  because 
of  v  dependence. 

Wait  and  Spies  (1964)  have  numerically  calculated  the  relationship 


between  changes  in  ;.l  and  h  and  changes  in  the  amplitude  and  phase  of 


VLF  signals  propagating  in  the  waveguide.  Ttieir  results  are  used  in 
linearized  expressions  for  computing  perturbations  of  the  D-region. 

Using  Taylor  series  expansion,  small  changes  in  amplitude  and 
phase  can  be  written  as 


As 


Ap 


=  Ah 


ll  * 

ah 

,  .  35 

A  t.' 

'  "  3f  ’ 

(2) 

+ 

ah 

AP  3-P- 
k  3P  ’ 

(3) 

where  s  is  signal  strength  in  dB,  p  is  signal  phase  in  radians,  and 
A  denotes  small  changes  about  mean  values  s,  p,  h,  or  £.  The  numerical 
values  for  the  coefficients  in  Equations  (2)  and  (3)  are  given  in  Table 
1 . 


Table  1 

Waveguide  Propagation  Parameters 


Parameter 

Day 

Night 

Average  of  c 
and  night 

3  (km  1  ) 

0.3 

0.54 

0.42 

h ( km  ^ ) 

70.0 

85.4 

77.7 

s(dB) 

-36.36 

-23.9 

30.13 

3S/3p(dB-km) 

58.1 

47.9 

53.0 

as/3h( dB/km) 

0.27 

-0.42 

-0.07 

p( radi ans ) 

5556.4 

5573.0 

5564.7 

3  p/ 3  p(  rad -kill) 

16.83 

16.91 

16.87 

3p/3h(rad/km) 

0.659 

0.660 

0.660 
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The  parameters  in  Table  1  were  determined  as  follows.  The  daytime 
values  of  and  h  were  assumed.  The  nighttime  values  were  computed  for 


the  average  day-to-night  variations  in  the  amplitude  and  phase  data. 

The  values  for  s,  p,  and  the  derivatives  were  obtained  from  Wait  and 
Spies  (1964)  assuming  propagation  of  only  the  first  (n  =  1)  waveguide 
mode.  Waveguide  attenuation,  phase  velocity,  and  mode  excitation  were 
considered  in  the  computations.  The  amplitude  of  the  second  (n  =  2) 
waveguide  mod°  is  at  least  21  dB  below  the  amplitude  of  the  first  mode 
for  the  14,250  km  propagation  path.  Crombie  (1958)  proposes  that,  for 
short  paths  greater  than  15,000  km>  long-path  propagation  may  become 
important.  In  our  case,  the  long  (25,770  km)  path  signal  amplitude  is 
calculated  to  be  at  least  18  dG  below  the  short  path  signal  amplitude. 
Consequently,  high  order  waveguide  modes  and  long-path  propagation  are 
neglected  in  this  work. 

Solving  for  Ah  and  A 6  in  Equations  (2)  and  (3)  yield  matrix  equations 
of  the  form 


Ah 

_J11 

1 

C\J 

AS 

.A|!. 

^21 

J^2 

Ap 

(4) 


The  values  for  the  matrix 


coefficients  are 


given  in  Table  2. 


6 


Table  2 


Coefficients  Relating  Amplitude  and  Phase 
Fluctuations  to  Changes  in  the  D-Region 


Day 

Night 

Average 

J11 

-0.50 

0.423 

0.466 

J12 

1  .72 

1.198 

1  . 466 

J21 

0.0193 

0.0165 

0.0182 

0.008 

0.0105 

0.0019 

LUNAR- INDUCED  VARIATIONS 

The  lunar-induced  semi-diurnal  (L^)  component  is  identified  by 
spectral  analysis  of  the  data.  A  time  series  representation  of  the 
component  is  provided  by  digital  filters.  The  filtering  techniques 
used  here  are  described  in  detail  by  Bernhardt  (1979). 

The  tidal  frequency  has  a  theoretical  value  of  1.932  cycles 
per  day.  Fourier  transforms  of  both  the  phase  and  amplitude  data  show 
a  spectral  peak  at  1.93  •  0.0035  cycles  per  day.  To  test  the  signifi¬ 
cance  of  that  frequency  line,  the  phase  data  were  divided  into  three, 
nonoverlapping  time  series  of  88  days  each.  The  Fourier  transform 
spectrum  of  each  individual  series  also  showed  the  line.  Figure  5 
illustrates  the  average  of  the  three  88  day  spectrums.  The  error  bar 
indicates  the  standard  deviation  of  the  L0  component. 

The  data  are  dominated  by  strong  solar  harmonics  having  frequencies 
that  are  multiples  of  1  cycle  per  day.  These  are  labeled  S-|  and  S., 
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in  Figure  5.  These  harmonics  are  suppressed  by  processing  the  data 
with  a  digital  filter.  Figure  6  illustrates  the  spectrum  of  the  VLF 


amplitude  data  after  the  filtering.  The  lunar  (L^)  spectral  line  is 
weak  in  comparison  with  other  spectral  lines.  Further  processing  is 
required  to  isolate  the  lunar  component. 

Both  the  amplitude  and  phase  data  are  next  filtered.  Only  fre¬ 
quencies  between  1.892  and  1.972  remain  in  the  data.  The  filter  used 
for  this  processing  has  a  response  time  of  25  days.  The  outputs  of 
the  filter  for  day  numbers  112  through  118  of  1967  are  illustrated  in 
Figure  7.  The  lunar-induced  phase  and  amplitude  are  described  by  the 
equations : 


APL  =  oip  cos 


2tt(x  -  6p) 

12 


radians  , 


(5) 


AS. 


ct_  cos 


2tt(t  -  Ss) 


12 


dB 


(6) 


where  x  is  the  local  lunar  hour,  x  =  0  hr  is  the  time  of  lower  lunar 
transit,  x  =  12  hr  is  the  time  of  upper  lunar  transit;  Ap^  is  the  lunar- 
induced  phase  variation;  As^  is  the  lunar  amplitude  variation;  ap  and  «s 
are  the  cosinusoidal  amplitudes;  and  6p  and  are  the  time  references 
in  lunar  hour.  The  lunar-induced  variations  have  the  seasonal  variations 
shown  in  Figure  8.  The  maximum  values  for  both  Ap^  and  As^  occur  near 
day  115.  The  reference  times  for  the  phase  and  amplitude  are  found  to 

be  &  =  5.48  ±  0.65  lunar  hour  and  6  =  8.42  +  0.16  lunar  hour. 

P  s 

The  amplitude  and  phase  changes  are  translated  into  changes  in 
the  D-region  by  using  Equation  (4).  Since  the  lunar  variations  are 


representative  of  both  the  day  and  night  D-region,  the  day-night 
average  values  for  J^,  an^  dj2  are  use^-  This  analysis 

yields  the  following  description  of  changes  in  the  D-region  due  to 
the  lunar  time  near  the  time  of  maximum  effect: 


AhL  =  0.39  cos[^-V^]  km  , 

(7) 

A3L  =  0.0048  cosp- (T  km'1  . 

(8) 

The  D-region  reference  height  reaches  a  maximum  value  0.39  km  above 
its  average  value,  at  approximately  5.9  and  17.9  lunar  hour.  The  D- 
region  bottom-side  scale  gradient  maximizes  near  8.2  and  20.2  lunar 
hour. 


HEBDOMADAL  VARIATIONS 

Spectral  analysis  of  the  phase  and  amplitude  data  does  not  produce 
a  clear  spectral  line  with  a  period  of  seven  days.  Thus,  hebdomadal 
effects  are  not  as  pronounced  as  the  lunar  tidal  effects.  This  may  be 
partially  due  to  the  increased  level  of  background  noise  at  the  lower 
frequencies. 

The  search  for  a  weekly  periodicity  in  the  VLF  phase  and  amplitude 
data  proceeded  as  follows.  The  data  were  separated  into  nighttime 
(2030-0350  local  time)  and  daytime  (0830-1550  local  time)  blocks.  The 
periods  were  chosen  to  exclude  the  transition  periods  at  dawn  and  dusk. 
The  data  were  then  superposed  and  averaged  according  to  the  day  of  the 
week.  Only  unbroken  strings  of  data  of  seven  days  were  used.  The 


superposition  process  is  equivalent  to  filtering  the  data  with  a 
harmonic  filter  which  passes  frequencies  which  are  multiples  of  1//’ 
cycle  per  day  (Bernhardt,  1974;  Bernhardt  and  Schlapp,  1976). 

The  results  of  the  superposition  analysis  are  illustrated  in 
Figure  8  and  Figure  9.  The  phase  does  not  show  any  significant  varia¬ 
tion  with  day  of  the  week.  At  night,  the  average  phase  on  Tuesdays 
and  Wednesdays  exceeds  the  average  phase  for  the  rest  of  the  days  of 
the  week  by  about  35  degrees.  This  is  not  significant  when  compared 
to  the  standard  deviation  of  the  measurements  which  is  greater  than 
70  degrees.  Similarly,  the  35  degree  drop  in  phase  during  the  Monday 
and  Tuesday  daylight  hours  is  not  considered  to  lie  significant. 

Analysis  of  the  amplitude  data,  however,  shows  a  clear  indication 
of  an  increase  in  VLF  amplitude  on  Sundays  (Figure  9).  At  night,  the 
amplitude  on  Sundays  is  0.14  dB  greater  than  the  weekly  average;  during 
the  day,  that  amplitude  is  0.46  dB  above  the  average.  These  increases 
are  almost  two  standard  deviations  above  the  mean. 

The  increase  in  amplitude  on  Sundays  cannot,  by  itself,  be  used 
to  determine  the  changes  in  the  h  and  6  of  the  D-region  model.  Using 
F.quation  (3)  and  Table  1,  one  obtains  the  expressions 


where  the 
night  and 


As^  =  0.42  Ah^  +  47.93  ,v^  =  0.14  dB  ,  Night 

ASp  =  0.27  .‘dip  +  58.1  .'.pg  =  0.46  dB  ,  Day 

subscripts  N  and  D  denote  hebdomadal  variations  during 
the  day,  respectively. 


(9) 

(10) 


the 
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DISCUSSION  or  RLSULTS 


The  analysis  presented  in  the  paper  has  shown  that  VLF  phase  and 
amplitude  measurements  show  evidence  of  lunar  tidal  oscillations.  The 
VLF  amplitude  data  also  contain  a  7-day  periodicity  which  may  be  related 
to  hebdomadal  variations  in  electric  power  consumption.  No  other  varia¬ 
tions,  such  as  might  be  produced  by  the  solar  rotation,  were  identified 
in  this  study. 

As  discussed  in  the  Introduction,  the  lunar-induced  variation  in 
VLF  phase  has  been  previously  observed  by  Brady  and  Crombie  (1963). 

They  concluded  that  the  D-region  reaches  a  peak  altitude  0.13  km  above 
its  average  altitude  at  a  lunar  hour  of  11.33  local  lunar  time.  The 
study  presented  here  gives  a  maximum  ,\hf  of  0.39  km  at  5.9  lunar  hour. 

The  factor  of  three  discrepancy  in  ,\h^  is  probably  due  to  seasonal 
variations  in  the  lunar  effect.  Some  differences  between  the  two  results 
may  also  be  due  to  the  different  paths  on  which  the  measurements  were 
made.  The  paths  used  by  Brady  and  Crombie  covered  more  than  26''  in 
longitude.  The  present  study  used  a  propagation  path  which  did  not 
vary  in  longitude  by  more  than  1.26'.  Since  the  magnitude  of  the  lunar 
tide  is  strongly  dependent  on  local  time,  one  expects  to  see  the  largest 
effect  for  a  VLF  path  which  stays  at  a  constant  longitude. 

The  time  reference  for  the  lunar  semi-diurnal  component  is  approxi¬ 
mately  b.5  lunar  hours  different  for  the  two  studies.  Since  the  lunar 
"semi-day"  is  12  hours  long,  the  two  results  are  nearly  1U0  degrees  out 
of  phase.  The  results  presented  here  show  that  the  D-region  reaches  its 


minimum  height  0.1  hours  before  lunar  transit.  This  is  consistent  with 
the  analysis  of  Appleton  and  Beynon  (1949)  showing  peak  2  MHz  absorption 
1.10  hours  before  transit,  and  it  is  also  in  agreement  with  the  0.15  MHz 
absorption  measurements  of  Mitra  (1955)  showing  peak  absorption  2.3  hours 
before  transit. 

As  expected,  the  lunar  tide  produces  VLF  amplitude  changes  along 
with  the  phase  changes.  The  largest  value  of  amplitude  variation  is 
.*0.26  dB.  The  amplitude  fluctuation  lags  the  phase  fluctuation  by  about 
three  hours.  A  *0.0048  km  ^  variation  in  the  normalized  bottom-side 
gradient  {,')  of  the  0-region  is  estimated.  An  increase  in  h  and  ,s 
probably  results  from  a  decrease  in  electron  concentration  and  an 
increase  in  the  neutral  concentration  at  the  bottom  of  the  ionosphere. 

The  results  of  our  analysis  can  be  interpreted  as  a  change  in  the  neutral 
concentration  due  to  the  lunar  semi-diurnal  tide. 

The  hebdomadal  variations  in  VLF  amplitude  are  consistent  with 
theories  that  attribute  precipitation  of  energetic  electrons  from  the 
magnetosphere  to  power  line  radiation  (PLR).  The  results  of  the  pre¬ 
ceding  section  show  that,  on  Sundays,  the  VLF  amplitude  is,  on  the 
average,  higher  than  during  the  rest  of  the  week.  This  is  most  likely 
due  to  some  anthropogenic  activity.  One  explanation  is  that  less  PLR 
occurs  on  Sunday  because  of  reduced  weekend  usage  of  electric  power  in 
the  western  United  States.  The  lower  level  of  PLR  in  the  magnetosphere 
scatters  less  radiation  belt  electrons  into  the  atmosphere.  This  would 
lead  to  a  decrease  in  the  D-region  ionization  resulting  in  an  increase 
in  the  parameters  h  and  r1.  According  to  Equations  (9)  and  (10),  this 
process  could  account  for  the  Sunday  increase  in  VLT  signal  strength. 
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FIGURE  CAPTIONS 


Figure  1.  VLF  propagation  path  between  NLK,  at  Jim  Creek,  Washington, 
and  Byrd  Station,  Antarctica. 

Figure  2.  Diurnal  variations  in  the  relative  amplitude  of  the  VLF 

signals  received  at  Byrd.  The  day  of  the  year  is  plotted  in  universal 
time.  Note  how  that  the  day-night  transitions  are  steeper  during  the 
period  near  Equinox. 

Figure  3.  Phase  measurements  made  at  Byrd.  Slow  trends  in  the  data 
are  due  to  frequency  differences  between  transmitter  and  receiver. 

Figure  4.  Diurnal  variations  in  the  relative  phase.  Note  the  rapid 
changes  at  dawn  and  dusk  during  the  period  near  Equinox. 

Figure  5.  Power  spectrum  of  the  phase  measurements.  The  solar  diurnal 
harmonics  (S^  and  S^)  dominate  the  spectrum.  The  magnitude  of  the 
lunar  (L^)  line  is  1/80  of  the  S-j  magnitude. 

Figure  6.  Power  spectrum  of  the  VLF  amplitude  measurements  with  sup¬ 
pressed  solar  harmonics. 

Figure  7.  Output  of  a  digital  filter  with  passband  at  1.932  ±  0.04 
cycles  per  day.  The  lower  and  upper  transit  of  the  moon  are  desig¬ 
nated  by  L  and  U,  respecti vely. 

Figure  8.  Seasonal  variations  in  the  amplitude  of  the  lunar-induced 
fluctuations  in  VLF  phase  and  signal  strength.  Both  observations 
have  maximum  values  near  day  115  of  1967. 

Figure  9.  VLF  phase  superposed  and  averaged  by  the  day  of  the  week. 

The  variations  are  a  fraction  of  one  standard  deviation  of  the  mea¬ 
surements  and  consequently  are  not  significant. 


Figure  10.  VLF  amplitude  superposed  and  averaged  by  day  of  the  week. 


A  significant  increase  in  signal  strength  occurs  on  Sunday  during  both 
night  and  day. 
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Figure  5.  Power  spectrum  of  the  phase  measurements.  The  solar  diurnal 
harmonics  (S]  and  S?)  dominate  the  spectrum.  The  magnitude  of  the  lunar 
(L2)  line  is  1/80  of  the  S]  magnitude. 
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Figure?  6.  Power  spectrum  of  the  VLF  amplitude  measurements  with  sup 
pressed  solar  harmonics. 
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Figure  7.  Output  of  a  digital  filter  with  passband  at  1.932  .t  0 
cycles  per  day.  The  lower  and  upper  transit  of  the  moon  are  d 
nated  by  L  and  U,  respectively. 
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Figure  8.  Seasonal  variations  in  the  amplitude  of  the  lunar-induced 
fluctuations  in  VLF  phase  and  signal  strength.  Both  observations 
have  maximum  values  near  day  115  of  1967. 
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